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Abstract 


Geometry is present in physics at many levels, most prominently in the theory of 
gravity. In these introductory lectures geometrical concepts like manifolds, geodesics, 
curvature and topology are introduced as a tool to describe and interpret physical 
phenomena in space-time, including the gravitational field itself. The focus is not only 
on possible static structures, e.g. black holes, but also on dynamical effects like waves 
and traveling domain walls. Although quantum gravity itself is outside the scope of 
these lectures, it is briefly discussed how geometry can be used even in that context 
to characterise configurations dominating the path-integral in various circumstances, 
including typical quantum processes like vacuum tunneling. Many of these geometrical 
methods can also be used in other branches of physics by linking the dynamics of a 
system with the geometry of its configuration space; some examples are mentioned. 
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Chapter 1 


Gravity and Geometry 


1.1 The gravitational force 


Gravity is the most universal force in nature. As far as we can tell from ob- 
servations and experiments every object, every particle in the universe attracts 
any other one by a force proportional to its mass. For slow moving bodies at 
large distances this is a central force, inversely proportional to the square of the 
distance. As the action is reciprocal, and since according to Newton action and 
reaction forces are equal in magnitude, the expression for the gravitational force 
between two objects of mass Mı and Mọ at a distance R is then determined to 
have the unique form 


M,M 
T 2, (1.1) 


The constant of proportionality, Newton’s constant of gravity, has dimensions 
of acceleration per unit of mass times an area. Therefore its numerical value 


F=G 


obviously depends on the choice of units. In the MKS system this is 


G = 6.67259(85) x 107! m? kg! s7?. (1.2) 


It is also possible, and sometimes convenient, to fix the unit of mass in such 
a way that Newtons constant has the numerical value G = 1. In the natural 
system of units, in which also the velocity of light and Planck’s constant are 
unity: c= ÅA = 1, this unit of mass is the Planck mass mp: 


mp = Vhe/G = 2.17671 x 10-8 kg 


= 1.221047 x 10'°GeV/c’. 


(1.3) 


Newton’s law of gravity (1.1) is valid for any two massive bodies, as long as they 
are far apart and do not move too fast with respect to one another. In particular, 
it describes the motions of celestial bodies like the moon circling the earth, or the 


planets orbiting the sun, as well as those of terrestrial objects like apples falling 
from a tree, or canon balls in free flight. Ever since Newton this unification of 
celestial and terrestial mechanics has continued to impress people and has had a 
tremendous impact on our view of the universe. It is the origin and basis for the 
belief in the general validity of physical laws independent of time and place. 


1.2 Fields 


Although as a force gravity is universal, Newton’s law (1.1) itself has only limited 
validity. Like Coulomb’s law for the electrostatic force between two fixed charges, 
Newton’s law holds strictly speaking only for static forces between bodies at rest. 
Moreover, and unlike the electric forces, there are modifications at smaller, finite 
distances which can be observed experimentally. 

For example, if the gravitational force would have a pure 1/R? dependence, 
the orbits of particles around a very heavy central body would be conic sections: 
ellipses, parabola’s or hyperbola’s, depending on the energy and angular mo- 
mentum, in accordance with Kepler’s laws. The observation of an excess in the 
precession of the perihelion of the orbit of Mercury around the sun by LeVerrier 
in 1845, and improved by Newcomb in 1882 [1], was one of the first clear indica- 
tions that this is actually not the case, and that the gravitational force is more 
complicated. 

The exact form of the gravitational forces exerted by moving bodies is a 
problem with many similarities to the analogous problem in electrodynamics. 
The understanding of electrodynamical phenomena greatly improved with the 
introduction of the concept of local field of force. This concept refers to the 
following characteristics of electrodynamical forces: 


e the influence of electric charges and currents, and of magnetic poles, extends 
throughout empty space; 


e the force on a standard test charge or test magnet at any given time depends 
on its location with respect to the source of the field and the relative state 
of motion; 


e changes in the sources of the fields give rise to changes in the force on test 
objects at a later time, depending on the distance; the speed of propagation 
of disturbances in empty space is finite. 


In the case of electrodynamics this speed turned out to be a universal constant, 
the speed of light c. One of the most striking consequences of these properties, 
which follow directly from the mathematical description of the fields as expressed 
by Maxwell’s equations, is the existence of electromagnetic waves, local variations 
in the fields which transport energy, momentum and angular momentum over 
large distances through empty space. Maxwell’s predictions were magnificently 


verified in the experiments of Hertz, confirming the reality and characteristics of 
electromagnetic waves [2]. From these results it became clear, for example, that 
light is just a special variety of such waves. Therefore optical and electromagnetic 
phenomena have a common origin and can be described in a single theoretical 
framework. 

The concept of field has found other applications in physics, for example in the 
phenomenological description of fluids and gases. Treating them as continua, one 
can describe the local pressure and temperature as fields, with local variations in 
these quantities propagating through the medium as sound waves or heat waves. 
As in electrodynamics, also in gravity the concept of field has taken a central 
place. What Maxwell achieved for electromagnetism, was accomplished by Ein- 
stein in the case of gravity [3]: to obtain a complete description of the forces 
in terms of space- and time-dependent fields. This field theory is known as the 
general theory of relativity. One of the most remarkable aspects of this theory 
is, that it provides an interpretation of gravitational phenomena in terms of the 
geometry of space-time. Many detailed presentations of the geometrical descrip- 
tion of space-time can be found in the literature, including a number of books 
included in the references at the end of these lecture notes [6]-[14]. A general 
outline of the structure of the theory is presented in the following sections. 


1.3 Geometrical interpretation of gravity 


The first step in describing a local field, both in electrodynamics and gravity, 
is to specify the potentials. The force on a test particle at a certain location 
is then computed from the gradients of these potentials. In the case of gravity 
these quantities have not only a dynamical interpretation, but also a geometrical 
one. Whereas the electromagnetic field is described in terms of one scalar and 
one vector potential, together forming the components of a four vector! A (x), 
the gravitational field is described by 10 potentials, which can be assembled in 
a symmetric four tensor gu (£). This tensor has a geometrical interpretation as 
the metric of space-time, determining invariant space-time intervals ds in terms 
of local coordinates x": 


ds? = gu (£) dz”dz”. (1.4) 


In the absence of gravity, space-time obeys the rules of Minkowski geometry well- 
known from special relativity. In cartesian coordinates (ct, x,y,z) the metric of 
Minkowski space is 


Nu = diag(—1,+1,+1,+1). (1.5) 


We use greek letters u,v, ..., taking values (0,1,2,3), to denote the space-time components 
of four vectors and tensors. 


Then the invariant ds? takes the form 


ds? = nu dr“dz” = ~ed? + dr? + dy? + dz’. (1.6) 


For a test particle, there is a simple physical interpretation of invariant space- 
time intervals measured along its world line, in terms of proper time measured 
by a clock in the local rest frame of the particle. As world lines are time-like, 
any interval along a world line satisifies ds? < 0, and the proper time interval dr 
measured by a clock at rest with respect to the particle is given by 


edr? = —ds’. (1.7) 


The local matrix inverse of the metric is denoted by g” (æ): 


a (a) gale) = 8, (1.8) 


where 6\ is the Kronecker delta, and the Einstein summation convention has 
been used. This convention implies that one should automatically perform a full 
summation over any repeated upper and lower indices, like v on the left-hand 
side of this equation. 

{From the gradients of the potentials, i.e. of the metric tensor, one constructs a 
quantity known as the connection, with components ry given by the expression: 


1 
Pa” = so" (Ong + Ov Gur = Onduv) . (1.9) 


At this point this expression may simply be considered a definition. It is a useful 
definition because this combination of gradients of the metric occurs very often 
in equations in general relativity. The deeper mathematical reason for this stems 
from the transformation properties of expressions for physical quantities under 
local coordinate transformations. However, it is not necessary to discuss these 
aspects at this point. 

A physical argument for the importance of the connection is, that it gives 
directly the gravitational acceleration of a test particle in a given coordinate 
system by the following prescription, due to Einstein: the world lines of point- 
like test particles in general relativity are time-like geodesics of space-time, the 
latter being considered as a continuous pseudo-euclidean space? with metric guv. 
Time-like geodesics are lines along which the total proper time between two fixed 
events fdr is an extremum; such world lines are solutions of the equation of 
motion 


?The term pseudo-euclidean here refers to the fact that the metric is not positive-definite, as 
it would be in a true euclidean space; space-time has a lorentzian signature (—, +, +, +), which 
implies that a general metric g,, has one negative (time-like) and three positive (space-like) 
eigenvalues. 


Cat u dz” dx* 
dr? àv dr dr 


A geometrical interpretation of the connection is therefore that it defines geodesics. 


= 0. (1.10) 


As concerns its dynamical interpretation, there is some similarity of eq.(1.10) with 
the Lorentz force in electrodynamics in that the proper acceleration depends on 
the four velocity t” (we use an overdot to denote a derivative w.r.t. proper time). 
However, whereas the Lorentz force is linear in the four velocity, the gravitational 
acceleration depends quadratically on the four velocity. The coefficients of this 
quadratic term are given by the connection, which therefore directly determines 
the force. Notice, that the gravitational acceleration does not involve the mass 
of the test particle. This is the mathematical content of the famous equivalence 
principle, which states that in a fixed gravitational field all bodies fall with the 
same acceleration, independent of their mass. 

Eqs.(1.10) and (1.9) express the force on a test particle in terms of a field 
(the connection), and the field in terms of potentials (the metric). This provides 
the basis for the field theory of gravity formulated by Einstein in his theory of 
general relativity. 

What is still lacking here of course is a set of equations determining the fields 
or potentials in terms of given sources, like Maxwell’s equations allow one to 
compute the electromagnetic fields in terms of charges and currents. General 
relativity provides such equations; they are known as the Einstein equations, 
and they determine the gravitational fields arising from given distributions of 
mass-energy and momentum. These distributions, the sources of the gravita- 
tional fields, are described by a symmetric tensor T u(x) with 10 independent 
components, equal to the number of independent components of the metric; Ty, 
is called the energy-momentum tensor. In order to give a concise formulation of 
Einstein’s equations, it is useful to first introduce one more geometrical concept, 
that of curvature. 


1.4 Curvature 


Curvature is a important concept, because it characterizes the space-time geom- 
etry in a way which is basically coordinate independent. In order to illustrate 
the idea, we first present a heuristic discussion. Then we turn to mathematical 
definitions, which allow us to write Einsteins equations in a compact way with a 
direct and elegant geometrical interpretation. 

Curvature is defined by the properties of surfaces; it is a measure of their 
deviation from a flat plane. Curvature can be expressed conveniently in terms of 
the change of a vector when transported parallel to itself around a loop enclosing 
a two-dimensional surface element. In a plane, displacement of a vector parallel 
to itself will always result in the same vector when you return to the starting 


point: there is absolute parallelism, but no curvature. In a curved surface this is 
no longer true. Consider the example of a sphere, a surface of constant curvature. 
If a vector is moved southward from the north pole parallel to itself, along a great 
circle, it ends up on the equator cutting it at a right angle. Now transport it 
along the equator over an arc of angular extension 0, always keeping at a right 
angle. Finally, move it back to the north pole along a great circle, keeping it 
parallel to the great circle and to itself. It then ends up at the north pole rotated 
by an angle 0 with respect to its original orientation. Thus parallel transport of a 
vector along a closed loop on a curved surface does not leave the vector invariant 
in general. 

In the example of the sphere, the rotation of the vector is proportional to the 
distance traveled along the equator, as measured by the angle 6. This is also 
proportional to the area enclosed by the loop consisting of the two sections of 
great circles through the north pole and the arc of the equator. For a sphere with 
radius R, the change dV in a vector V when the loop encloses an area dA = R?d6 
is 


V - 
dV =Vdb = 73 dA = KV dA. (1.11) 


The constant of proportionality K = 1/R? between the relative change in the 
vector dV/V and the change in the area dA is a measure of the intrinsic curva- 
ture of the surface, known as the Gaussian curvature. Note, that the Gaussian 
curvature is large when the radius of the sphere is small, and that it vanishes in 
the limit R — oo, as one would expect intuitively. 


This idea can now be generalized to situations in higher-dimensional spaces. Of 
course, in a space of three or more dimensions, there are many surfaces one can 
draw through a given point. Each of these surfaces may have its own curva- 
ture. Our first conclusion is therefore, that in general the curvature has many 
components. 

Fortunately, it is only necessary to consider curvature components in surfaces 
which are linearly independent. To this end, let us choose a locally non-degenerate 
coordinate system x“ and consider in particular the surfaces defined by keeping 
all coordinates constant except two, say £” and x”. 

Let d£“ = —d&”" be the area of an oriented surface element in the x“-x”- 
plane which contains the point with coordinates (£o). Any vector with compo- 
nents Vò transported parallel to itself around the loop defined by the circumfer- 
ence of this surface element now changes to first order by an amount 


1 
2 
The constants of proportionality Rus (Xo) define a tensor known as the Riemann 
curvature tensor at the location (xo). A comparison with eq.(1.11) shows, that the 


dV* = —= R uk (£0) Ved". (1.12) 


Riemann tensor is a higher-dimensional generalization of the Gaussian curvature. 
From eq.(1.12) it follows immediately, that a space-time is completely flat if and 
only if all components of the Riemann curvature tensor vanish. 

Defining parallel transport of a vector along a line element dx* in terms of 
the Riemann-Christoffel connection (1.9): 


dV" = —dx* D" V”, (1.13) 
a direct computation for an infinitesimal closed loop gives the result that the 
components of the Riemann tensor can be expressed as 


Rae = Wig — OP ye — Pe EA, (1.14) 


u 


where the square brackets denote an ordinary matrix commutator of two connec- 
tion components, considered as 4 x 4 matrices (Ta), À, the first one with a = p, 
the second one with a = v. 

Using eq.(1.9) for the connection, we can ultimately express the curvature in 
terms of the metric and its first and second derivatives. It may then be checked 
from this explicit expression that it identically satisfies the relation 


Ruy” + Roya + Ryu = 0. (1.15) 


This relation is known as the Bianchi identity. It is analogous to the homogeneous 
Maxwell-equation, which implies that the electric and magnetic fields can be 
obtained from a four-vector potential A,. In a similar way the Bianchi identity 
for the curvature tensor implies that the curvature tensor can be obtained from 
a symmetric tensor guv- 

We wish to emphasize here, that the curvature of space-time is not a measure 
for the strength of the gravitational fields or forces; rather, it is a measure for 
the variation of the gravitational fields in space and time, the gradients. 


1.5 The Einstein equations 


The Einstein equations are second order non-linear partial differential equations 
for the metric of space-time in terms of a given energy-momentum distribution. 
As a result, these equations describe variations of the fields in space and time, 
rather than the fields themselves. In view of the role of curvature as a measure for 
such variations and the geometrical interpetation of gravitational forces, it is quite 
natural that the Einstein equations should admit an interpretation as equations 
for the curvature of space-time. However, rather than the full Riemann tensor, it 
turns out that these equations involve only its contracted forms, the Ricci tensor 
Ra and the curvature scalar R: 


Ry = Ry, R= RE = g Ry. (1.16) 


prAvY 9 


Note that we adhere to the general convention that contraction with an inverse 
metric can be used to raise an index of a vector or tensor, turning a covariant 
component into a contravariant one; similarly contraction of an upper index on 
a vector or tensor with the metric itself is used to lower it to make a covariant 
component out of a contravariant one. 

The reason that only the contracted forms of the Riemann curvature appear 
in the Einstein equations is, that in four-dimensional space-time the Riemann 
curvature tensor has too many components, twice as many as the metric. There- 
fore fixing the metric components in terms of the full Riemann tensor would 
generally lead to an overdetermined system of equations. In contrast, the Ricci 
tensor is symmetric: R,, = Rvp, and has 10 independent components, precisely 
equal to the number of metric components to be solved for. It is therefore to be 
expected that the Einstein equations take the form of an equation for the Ricci 
tensor. Indeed, Einstein’s equations for the gravitational fields can be written in 
the simple form 


1 8rG 
Ryu — zIw R = ZT tw» (1.17) 
or equivalently: 
8SrG 1 
Rw = -2 (Ta — 5 dw ) (1.18) 


where T = T,“ = g'”T,,. Notice the appearance of Newton’s constant on the 
right-hand side of the equations. Of course, in empty space-time Tp» = 0; the 
second version of Einstein’s equations (1.18) then implies that the Ricci tensor 
vanishes: 


Ry = 0. (1.19) 


Space-time geometries which satisfy these conditions are some-times called Ricci- 
flat. The most interesting aspect of these equations is, that an empty space-time 
can still support non-trivial gravitational fields: the Riemann curvature tensor 
can have non-zero components even if its contracted form, the Ricci tensor, is 
zero everywhere. 

This observation is important, because it lies at the heart of the theoretical 
arguments for the existence of gravitational waves. Indeed, as we discuss in 
more detail later on, one can find non-trivial solutions of the Einstein equations 
representing lumps in the gravitational field which travel at the speed of light 
and transport a finite amount of energy and momentum per unit of volume from 
one flat part of empty space to another flat part of empty space. 

There is a modified version of Einstein’s equations in which the gravitational 
field itself creates a non-vanishing energy-momentum tensor, even in the absence 
of matter. This gravitational field behaves like a homogeneous continuum with 
constant density and pressure (one of them negative), satisfying the equations 


8rG 
Ryu = Tr Iw, (1.20) 


with A a constant. Such an equation arises if the energy-momentum tensor is 


Tw = JuA. (1.21) 


The constant A is known as the cosmological constant. Contraction of eq.(1.20) 
with the inverse metric gives 


R = TG A (1.22) 


ct 


This shows that in the presence of a cosmological constant empty space-time has 
a constant non-zero scalar curvature. Depending on the sign of the cosmological 
constant, eq.(1.20) then admits solutions of Einstein’s equations corresponding 
to empty space-times in continuous expansion or contraction. This modified 
version of Einstein’s equations is especially interesting with regard to cosmological 
applications. 

Of course in the general case there is a contribution to the energy-momentum 
tensor from all matter and energy which happens to be present in the part of 
the universe we wish to describe. The precise form of the energy-momentum 
tensor depends on the kind of matter or radiation which constitutes the source 
of the gravitational field. However, it is very often justified and convenient to 
treat matter and radiation as an ideal fluid, the energy-momentum tensor of 
which depends only on the mass- and pressure-density and its flow, as decribed 
by the local four velocity of an element of fluid at location with coordinates (x). 
Denoting the mass density by p(x), the pressure density by p(x), and the local 
four velocity by u(x), the energy-momentum tensor then takes the form 


Tuv = PIu + (p+ p) upu. (1.23) 


Comparison with eq.(1.21) shows, that the cosmological constant can be reob- 
tained from an ideal fluid with p = —p = A. 


1.6 The action principle 


The Einstein equations (1.17) can be obtained from a variational principle; that 
is, they follow by requiring an action S to be stationary under variations dg, of 
the metric: 6S = 0. The action for the gravitational field equations, including a 
coupling to material energy-momentum, is 


ct 


— 1 4 — uv 
S = z [tev a (SRo o). (1.24) 


10 


Here O, (x£) represents the specific source term, but generally it is not the energy- 
momentum tensor itself; rather, the energy-momentum tensor is given in terms 
of ©,, and its covariant trace O = O,“ by 


Tw = Ow — 5 guð. (1.25) 
However, we consider 0, to be the independent variable, hence by fiat there 
is no contribution from this tensor under variations with respect to the metric. 
Moreover, in eq.(1.24) g(x) = det gu (x), and the minus sign is necessary because 
of the lorentzian signature of the metric, which implies the product of its eigen- 
values to be negative everywhere. The proof of the action principle can be stated 
in a few lines, but it requires some intermediate results which we discuss first. 
To begin with, in stead of considering variations of the metric g,,, one can 
equivalently use the variations of the inverse metric g"”. Indeed, because g,,,g°” = 
ôA, it follows that 


Sga g” + gag” = 0, (1.26) 


or 


O Guy = 4.59 Gav (1.27) 


Thus the variation of the metric can be expressed in terms of the variation of the 
inverse metric and vice versa. 

Next, we need a rule for computing the variation of the determinant g. This 
is obtained from the formula 


log det M = Tr log M, (1.28) 


which holds for symmetric matrices M, as can easily be checked by diagonalizing 
M and using the rule that the logarithm of the product equals the sum of the 
logarithms of the eigenvalues. From this equation we then derive 


1 
z3 $9) = g” guv = —6 g” guv. (1.29) 


Finally, we observe that the Ricci tensor R,» is an expression involving only the 
connection PÀ. Hence 


SR = bg!” Ruy + g” Ru, (1.30) 


with the variation ó R„, expressible purely terms of the variation of the connection 
ôT aò. The total variation of the action now becomes 


11 


(1.31) 
The last term is to be interpreted in the sense that we vary R with respect to 
the metric only via its dependence on the connection. However, it turns out that 
if we insert the expression (1.9) for the connection, this whole variation adds up 
to a total derivative. As the action principle applies to variations which vanish 
on the boundary, this term does not contribute to 6S. If we then require the 
variation of the action to vanish, we indeed obtain the Einstein equations: 


4 


c 1 1 
spg Pe gwk) + Ow — gow 


(1.32) 


ct 


1 


The equation with the cosmological constant is obtained by making the replace- 
ment g” O, = O —> 2A in the action. 


| 
> 
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Chapter 2 


Geodesics 


2.1 Curves and geodesics 


In this chapter we discuss how basic information about the geometry of spaces 
or space-times —generically manifolds— can be obtained from studying their 
geodesics. The geodesics of euclidean differentiable manifolds are smooth curves 
of stationary proper length, meaning that to first order the proper length of an 
arc does not change under small variations of the curve. More generally and more 
precisely, given a metric on a manifold (euclidean or lorentzian) the geodesics are 
those smooth curves for which the proper interval between any two fixed points 
P and P’ through which the curve passes is both well-defined and stationary. 
The interval among the fixed points can be either a maximum, a minimum or a 
saddle point. 

A direct reason for geodesics to be important is that time-like geodesics rep- 
resent the allowed trajectories of test particles in a fixed background space-time. 
But in general geodesics tell us interesting things about the symmetries and 
topology of a manifold, which often have wider implications for physics. 

First we derive the fundamental equation characterizing geodesics. The proper 
interval As [Cpp:] of the arc Cpp of an arbitrary smooth curve C between the two 
fixed points (P, P’) is given by the integral 


da” da” 
Iw Iy dA) 


(2.1) 


As [Cpp] = [ _ t= [ 


This holds for all geodesics of euclidean manifolds, and time- or space-like geodesics 
of lorentzian manifolds; À is a continuous parameter on some interval of the real 
line, parametrizing the curve between the points P and P’ in any convenient way, 
provided the coordinates z”(À) of points on the curve are monotonic differentiable 
functions of this parameter: the label of any point must be unique and the tan- 
gent vector t“ = dx" /dX must be well-defined. The parameter À is then called an 
affine parameter. Note that the integral does not depend on the particular choice 


PP! 


13 


14 


of affine parameter: it is invariant under differentiable changes of parametriza- 
tion A + A’‘(A). This freedom to choose for A any smooth parametrization of 
the curve is often useful, but requires some care in interpreting results explicitly 
containing A. As any invariant geometrical results with physical interpretation 
can not depend on the choice of affine parameter, the best way to proceed is to 
remove any A-dependence from equations before drawing physical conclusions. 

The proper interval A s[Cpp/] is stationary if the integral does not change to 
first order under arbitrary smooth variations of the path while keeping the end 
points of the integral fixed: 


ô dà wļe(à)] = 0, (2.2) 
Cpp! 
where the integrand is w[a(A)] = 4/|gut”z"|. Now after a partial integration 


with vanishing boundary terms (the end points remain fixed) the variation of the 
integrand can be written 


d | guv dx” 1 gy dx* dz” 
7 { — hy — —_ 2.3 
o (5 [eee] + pee (2.3) 


provided w # 0. The sign on the r.h.s. is positive in euclidean spaces or for 
space-like curves in lorentzian space-times, and negative for time-like curves in 
lorentzian space-times. Light-like curves in lorentzian space-time need special 


ow = 


consideration, as ds? = 0 is equivalent to wdà = 0. 

The variational principle (2.2) now implies that for any smooth local varia- 
tion da" the induced variation dw must vanish if the curve «“(A) is to describe 
a geodesic. For w 4 0 we can divide the right-hand side of eq.(2.3) by w once 
more, and use the fact that wd\ = ds, the infinitesimal proper interval which 
is measured in the laboratory frame and independent of the choice of affine pa- 
rameter A. We can also replace the partial derivative of the metric in the second 
term by the connection; recall, that the connection was defined in eq.(1.9) as 


l gyw Irv gyv 
Ta = 5g" - . 2.4 
wv 24 ( Ox” Ox ðr“ (2.4) 
This expression is actually the solution of the equation 
Dgw 89 ; . 
Deo = Aw ~~ Py Gav ~~ ry, Jur = 0, (2.5) 


stating that the metric is covariantly constant, a condition known as the metric 
postulate. From the equation above it follows that 


gyv 
ðr” 


= Laa gew + Diy gra. (2.6) 


15 


The vanishing of dw is now seen to imply a differential equation of the form (1.10) 
for the coordinates of a point moving on the curve as a function of the proper 
interval parameter s: 


Cat dx’ dz" 
a4 u 
ds? *’ ds ds 


Note, that the above procedure eliminating the affine parameter 4 is equivalent 


= 0. (2.7) 


to a choice of rate of flow of the affine parameter dà equal to the proper interval 
ds measured in the coordinate frame {x}. If the signature of the manifold is 
lorentzian and the curve time-like this affine parameter is the proper time T, and 
eq.(2.7) is the equation of motion of a test particle of arbitrary non-vanishing 
mass, in a fixed space-time with affine connection [,“(a). Note also that the 
choice of dà = ds implies an identity 


dx" dx’ 
Tuy ds ds 


where again the plus sign holds for euclidean or space-like lorentzian curves, and 


= +1, (2.8) 


the minus sign for time-like curves in lorentzian space-times. This is a constraint 
on the solutions of the geodesic equation (2.7). 

That this constraint is consistent with the geodesic equation follows by ob- 
serving that 


1 dz” dx” 
H(s) = = gu ——— 2. 
(s) g Iu ds ds ( 9) 
is a constant of motion, or more generally a constant of geodesic flow: 
dH | dx" [dz 4 „dzò dz" _ 9 (2.10) 
ds ds \ ds “ds ds} ` 
Finally light-like curves (also called null-curves) by definition obey 
dz” dx” 
„— = = 0. 2.11 


In this case the affine parameter À can not be identified with proper time or 
distance. As eq.(2.11) is to hold on all points of the light-like curve, the quantity 


o= 1l dat dx’ 
~ ZIA aX 
must be a constant of light-like motion. Therefore eq.(2.10) still holds with s 
replaced by à. From this it follows that it is correct to interpret the solutions of 
the geodesic equation 


H(A) 


=0 (2.12) 


Cat dx” da" _ 9 (2.13) 
dd? weds dy” ` 
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satisfying the constraint H(A) = 0 as light-like geodesics. 


Remark. Although we have shown that the geodesic equation implies the con- 
stancy of H(A) for all types of geodesics, the inverse statement is not true. In par- 
ticular, for non-null geodesics (euclidean, and space-like or time-like lorentzian) 
the more general condition 


Pa” dx” dz" dx” 

— K = B* (zr)— 2.14 

I t ay dn Ae) (2.14) 
with B,, = —B,, an anti-symmetric tensor, is sufficient for H to be conserved. 


This happens for example for charged particles in an electro-magnetic field with 


field strength Fux: 


q 
Bywl(«) = = Ful), (2.15) 
m 
where q/m is the charge-to-mass ratio of the particle. Similar extensions exists 
for particles with spin or non-abelian charges. 
For light-like geodesics the equation can be further generalized to read 


Cx" dz” dx" dz” dz” 
ee pie EO Bl (gp) oe 
De + Gy aA vay + ka 


with x an arbitrary constant. It should be noted that if x 4 0, then H is only 
conserved by virtue of the fact that it vanishes. 


(2.16) 


Summary 
The results of this discussion an be summarized as follows. All geodesic curves, 
in euclidean or lorentzian manifolds, satisfy the equation 


Par” dx” dz" 
2 H = 2.1 
ae tt Gyan (2-17) 
where À is an affine parameter. The quantity 
1 dz” dx” 
HOA) = 5 Iw Gy (2.18) 


is conserved by the geodesic flow. In physical applications three cases are to be 
distinguished: 


1. For euclidean manifolds and space-like geodesics of lorentzian manifolds one 
can take dà = ds, the proper distance; these geodesics are then character- 
ized by the condition 2H(s) = 1. 


2. Time-like geodesics of lorentzian manifolds can be parametrized by the 
proper time: dà = dr; such geodesics satisfy 2H(7) = —1. 
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3. Light-like (null) geodesics admit no such identification of A with invariant 
measures of the arc length in an inertial frame (‘laboratory frame’); null- 
geodesics are characterized by H(A) = 0. 


2.2 Canonical formulation 


In D-dimensional space-time eq.(2.17) represents a set of D second order differ- 
ential equations in the affine parameter À providing a covariant description of 
geodesics. It is possible to replace this set by a set of 2D first order differential 
equations by going over to the canonical formulation of the dynamics of test par- 
ticles, in which the number of independent variables is doubled. Introducing the 
momentum variables 


dx” 
Palà) = Iw (2); (2.19) 
the geodesic equation (2.17) becomes 
dp z 
P = TAP. (2.20) 


The equations (2.19) and (2.20) constitute a pair of first-order differential equa- 
tions equivalent to the second order geodesic equation (2.17). A powerful result 
is now obtained by observing that these two equations can be derived from the 
constant of motion H(A), eq.(2.18) as a pair of canonical hamiltonian equations 
with H(A) in the role of hamiltonian: 


1 
H(A) = 5 9" PuPv- (2.21) 
With this definition 
u 
de" oH dpa _ _ OH (2.22) 
dà opu dà Ox" 


Then any function F(x", p,) on the 2D-dimensional phase space spanned by the 
coordinates and momenta changes along a geodesic according to 


dF dx" OF dp, OF 
a Pw OF _ ep 2.2 
Tx ~ dh Oo + WD Op, > OMI (2.23) 


where the Poisson bracket of two functions (F, G) on the phase space is defined 
generally as 


OF 0G OF OG 
Ox" Op, Öp, Ba" 
The anti-symmetry of the Poisson bracket automatically guarantees the conser- 
vation of the hamiltonian H by 


{F,G} = —{G,F} = (2.24) 
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{H,H} = 0. (2.25) 
A similar property holds for all quantities which are constant along geodesics: 
dF 
dà 
Formal properties of the Poisson bracket include its anti-symmetry and its lin- 
earity in each of the arguments: 


0 7 s {F H}=0. (2.26) 


{F aG + aG} = Q1 {F, Gi} +a {F, Ga}. (2.27) 
Another important property is that it satisfies the Jacobi identity: 


{F {G,K} +1G,{1K, FH} +{K, {F,G}} = 0. (2.28) 


This property implies, that the Poisson bracket of two constants of geodesic flow 


is again a constant of geodesic flow: if {F, H} = 0 and {G, H} = 0, then 


{{F, G} ’ A} = 0. (2.29) 


As H is itself constant on geodesics, the above results are sufficient to establish 
that the constants of geodesic flow form a Lie algebra, with the Poisson bracket 
as the Lie bracket. 

The definition of the Poisson bracket (2.24) is not manifestly covariant, and 
therefore at first sight it seems to destroy the build-in covariance of the tensor 
calculus of differential geometry, and of general relativity in particular. However, 
note that in our definition the connection (2.4) is alway symmetric in its lower 
indices; therefore we can rewrite the Poisson bracket formula as 


OG OF 
FG} =D,F— —- —D,G, 2.30 
{ } H p, p, H ( ) 
where for scalar functions F and G on the phase space we define the covariant 
derivative 


OF 
DF = ð F +1’ p—. (2.31) 
Op, 
This equation preserves manifest covariance for scalar functions F and G; for 
example, if J(x, p) = J"(x)p,, then 


Dit = (,J” + PJ) pv = (DS) pv. (2.32) 


The manifestly covariant form of the Poisson bracket can also be extended to all 
completely symmetric tensors T,,...,,, of rank n by contracting all n indices with 
the momentum p, to obtain the scalar 
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1 
T(x, p) = zy T (x) Pu Pans (2.33) 


Inserting this scalar into the covariant Poisson bracket, the bracket for a particular 
tensor component 7#1-#» is obtained by taking, in the expression resulting on 
the right-hand side, the coefficient of the n-nomial in the momenta that one is 
interested in. For anti-symmetric index structures a similar result can be achieved 
using differential forms, or more generally Grassmann algebras [18, 19, 20]. 


2.3 Action principles 


The canonical formulation of geodesic flow presented here can be cast into the 
form of a standard hamiltonian action principle. Namely, the canonical eqs.(2.22) 
define the critical points of the phase-space action 


2 
S= f dà (è - p — H(z, p)). (2.34) 

1 

By using the hamilton equation 


pt — OH 

OD 
the momentum variable can be eliminated in favour of the velocity of geodesic 
flow, to give the lagrangian form of the action 


2 | 
s=| ANS Juri a)dhe”. (2.35) 
1 


We observe, that this is not the action for the geodesics that we started from, 
eq.(2.1). In particular, it is not reparametrization invariant. A comparison shows, 
that the lagrangian in (2.35) is actually the square of w(A) used earlier, eq.(2.2). 
Clearly, the critical points are therefore the same, with the possible exception of 
light-like curves. Actually, the variational principle for the original action was 
not even well-defined in the case of light-like curves. 

In this section we examine somewhat more closely the relation between the 
various action principles, and in particular the role of reparametrization invari- 
ance. It makes clear how one can derive the lagrangian and hamiltonian actions 
presented above within a unified framework. 

The key to the general formulation is to incorporate reparametrization invari- 
ance in an action of the type (2.35) [21]. This is achieved by introduction of a 
new variable e(A), a function of the affine parameter, which acts as a gauge field 
in one dimension for local reparametrizations of the worldline. More precisely, 
the variable e(A) by definition transforms under reparametrizations A — A’ as 
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a 
dA 


This is the transformation rule for the square root of a metric in one dimension 


e(A) 4 e'(X') = e(A) (2.36) 


(on a curve), for which reason the variable e(A) is called the einbein. Written in a 
slightly different way, eq.(2.36) states that e(A)dA is defined to be invariant under 
reparametrizations. This makes it straightforward to write a reparametrization- 
invariant form of the quadratic action (2.35): 


2 


2 1 dz“ dz” œ 
See] =f dd (Same = -Ze (2.37) 


The last term, proportional to the constant c°, does not depend on the coordinates 


x". However, it is reparametrization invariant and its inclusion is important for 
reasons given below. 

Variation of the full action Sle; 2] w.r.t. the coordinates x" reproduces the 
geodesic equation (2.17) provided everywhere the differentials d\ are replaced by 
edX. This replacement makes the geodesic equation reparametrization invariant, 
as ed\ = e'dà'. On the other hand, the substitution e(A) = 1, which breaks 
reparametrization invariance, makes the old and new equations fully identical. 
Such a substitution is allowed because it only serves to define the affine param- 
eter A in terms of quantities measurable by an observer in the locally euclidean 
or lorentzian laboratory frame, like proper length or proper time. Of course, 
such a choice necessarily breaks the reparametrization invariance, but it pre- 
serves the physical content of the equations. It amounts to a choice of gauge 
for reparametrizations in the true sense of the expression, with the result that 
modulo a constant the action S[e; x] reduces to the simple quadratic action (2.35). 

However, the action $[e; 2] contains more information, as one can also consider 
its variation w.r.t. the einbein e. Requiring stationarity under this variation gives 
a constraint on the solutions of the geodesic equation, which after multiplication 
by 2e? takes the form 


Jul LJE E” + ee = 0. (2.38) 


With e = 1, comparison with the definition of H in eq.(2.9) shows that this 
constraint expresses the constancy of the hamiltonian along a geodesic: 


2H(A) = —e. (2.39) 
For euclidean manifolds one can always choose units such that ¢? = —1, but for 
lorentzian space-times c? = (—1,0, +1) corresponds to space-like, light-like or 


time-like geodesics respectively; in the latter case the constant c then represents 
the speed of light, which is unity when expressed in natural units. 

Finally we derive the original reparametrization-invariant geodesic action (2.1) 
from Sle; z]. This is achieved by using eq.(2.38) to solve for e, instead of fixing a 


21 


gauge: 


e = +y/ [gu thi”/c?|. (2.40) 


Substitution of this expression for e into the action Sle; z] then leads back (mod- 
ulo a possible sign) to the original action (2.1). Therefore we conclude that both 
this action and the quadratic action (2.35) follow from the same starting point 


provided by eq.(2.37). 


2.4 Symmetries and Killing vectors 


Often space-times of interest possess special symmetries. Such symmetries are 
useful in finding solutions of the geodesic equations, especially if they are contin- 
uous symmetries as these generate constants of geodesic flow. This is the content 
of Noether’s theorem applied to the motion of test particles and its generalization 
to space-like curves; see for example refs.[11]-[15], [18, 19]. In this section the role 
of continuous symmetries is explored in some detail. 

A manifold possesses a continuous symmetry, or isometry, if there are in- 
finitesimal coordinate transformations leaving the metric invariant. Finite trans- 
formations can be obtained by reiteration of the infinitesimal transformations, 
and donot require separate discussion at this point. Consider an infinitesimal 
coordinate transformation characterised by local parameters é“ (æ): 


rh — g" = g" — EK (e). (2.41) 


For the line element ds? to be invariant the metric must transform as a tensor; 
therefore by definition 


, ; ðr!" Ox" 
JuhE) = Ina (2.42) 


On the other hand, invariance of the metric under a coordinate transformation 
implies 


Gu (©) = Ju (2). (2.43) 
Combining the two equations (2.42) and (2.43) to eliminate the transformed 


metric g, and inserting the infinitesimal transformation (2.41), expansion to 
first order in €" leads to the covariant condition 


Dye, + DLE, = 0. (2.44) 


This equation is known as the Killing equation, and the vector fields defined by 
the solutions €"(2) are called Killing vectors. The existence of a Killing vector is 
equivalent with the existence of an infinitesimal coordinate transformation (2.41) 
leaving the metric invariant. 
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Let us pause for a moment to reflect on the interpretation of this result. If the 
coordinate transformation is considered to be a passive transformation, the new 
coordinates represent a relabeling of the points of the manifold: x“ and x’ in 
eq.(2.41) refer to the same point P. Now consider the point P’, which in the new 
coordinate system has the same value of the coordinates as the point P in the 
old coordinate system: «/"(P’) = «"(P). Eq.(2.43) states that the new metric at 
the point P’ has the same functional dependence on the new coordinates as the 
old metric on the old coordinates at the point P. Thus in the new coordinate 
system the neighborhood of P’ looks identical to the neighborhood of P in the old 
coordinates. In particular it is possible to map the points in the neighborhood of 
P one-to-one to points in the neighborhood of P’ in such a way that all distances 
(intervals) are preserved. As this is in fact a coordinate-invariant statement, 
it implies a local indistinguishability of the two neighborhoods of the manifold 
(as long as it is not endowed with additional structures). This then is what is 
implied by a continuous symmetry of the manifold generated by the solutions of 
the Killing equation. 

Now if two neighborhoods of the manifold are similar in the sense described 
above, then this similarity should also extend to the geodesics passing through 
the points mapped on each other by the symmetry operation. Indeed, as the new 
metric at the point P’ has the same functional dependence on the new coordinates 
there as the old metric at P has on the old coordinates, and not just the same 
value, the connections and curvature components there also become identical after 
the coordinate transformation. The same holds for the solutions of the geodesic 
equations. Thus the symmetry transformation maps geodesics through the point 
P to geodesics through the point P’, and the geodesic flow in the neighborhood 
of these points is completely similar. 


The similarity of geodesic structure generated by symmetry operations implies 
the conservation of certain dynamical quantities under geodesic flow. This is most 
conveniently discussed in the canonical formalism of the geodesic flow presented 
above in sect.2.2. The quantity conserved is the generator of the symmetry 
transformations, obtained by contraction of the Killing vector with the canonical 
momentum: 


JIE] = O(a) Dy (2.45) 


This is a scalar quantity, hence its value is coordinate invariant. It generates the 
infinitesimal coordinate transformations through the Poisson bracket (2.24): 


bu" = —& (x) = {J[e],«"}. (2.46) 


One can similarly determine the variation of the momentum under the transfor- 
mation generated by J[€]: 
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OP, = {JE Pub 


= oJ] = T,"(x)pv 


with Z,” = 0,€”. This represents a coordinate-dependent linear transformation 


(2.47) 


of the momentum. 

That the quantity J[é| is conserved now follows from the observation that 
the hamiltonian is invariant under the phase-space transformations defined in 
eqs.(2.46), (2.47). Indeed, taking into account the canonical equations of motion 
(2.23) and (2.26), and using the Killing equation (2.44) we obtain 


aig) _ E 
OM = {JH} = 0. (2.48) 


This result was to be expected, as the hamiltonian depends on the coordinates 
only through the metric, which is invariant by definition of an isometry. 


As the Killing equation (2.44) is linear it follows that the Killing vectors define 
a linear vector space: any linear combination of two Killing vectors is again a 
Killing vector. Let the dimension of this vector space be r; then any Killing 
vector can be expanded in terms of r linearly independent basis vectors e;(x), 
i= 1l,... r: 


Elai) = ayer(x) +... + ae, (x). (2.49) 


Equivalently, any infinitesimal symmetry transformation (2.41) depends linearly 
on r parameters a;, to which correspond an equal number of conserved generators 


JIE] = ay dy +... + oy JS), where J;(a,p) = ef (@)py. (2.50) 


We now show that these generators define a Lie algebra through their Poisson 
brackets. The key observation is, that the Poisson bracket of any two constants 
of geodesic flow is itself a constant of geodesic flow. This is implied by the Jacobi 
identity for the Poisson bracket, which can be written as: 


{{F,G},H} = {{F,H},G} — {{G, H}, F}. (2.51) 


Thus, if the Poisson brackets of F and G with the hamiltonian H vanish, then 
the left-hand side vanishes as well, and {F, G} is itself a constant of geodesic flow. 

A second important observation is, that for two constants of motion which 
are linear in the momentum, like the J;, the Poisson bracket is also linear in the 
momentum: 


La J} = (eð e! — dye") Pu- (2.52) 
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As the {J;}'_, are supposed to form a complete set, the expression on the right- 
hand side can be expanded in terms of the basis elements, with the result that 


with a certain set of coefficients f;; k> i we have 
{Ji, Ji} = fj" Jr. (2.53) 
An equivalent statement is that the basis vectors of the Killing space satisfy 
eOe; — ef 0e = oj eke (2.54) 


Eq.(2.53), combined with the observations that the J; form a linear vector space 
and that the Poisson bracket is bilinear and anti-symmetric in its arguments 
complete the proof that the generators of the symmetry transformations define a 
Lie algebra with structure constants f,,*. 

Of course the linear p.d.e.(2.54) encodes the same information. It can be 
interpreted in terms of the Lie-derivative, which is the operator comparing the 
value of a function at the point P with that at the point P’ which has the same 
value of its coordinates after an arbitrary transformation of the form (2.41) as 
the point P had before the transformation: «/"(P’) = «"(P). For a scalar ¢(x), 
for which ¢'(#’) = g(x), this gives 


Le [p] (x) = ling [Ø (x) — O(2)] = Eagle), (2.55) 


where the ling[Q] denotes the part of the expression Q linear in €. For a vector 


v(x), for which 


(2.56) 


one finds similarly 


(Le [v])" (2) = ling [v!(x) = v"(w)] = €°,v"(w) — déo (a). (2.57) 


For higher-rank tensors the construction of the Lie-derivative gives the result 


(Le [TJU (e) = ling [Pte e) — TH!) 


EO TEHE (x) — O E TE (x) = 1. = OET” (x), 

(2.58) 
Now comparing eqs.(2.54) and (2.57) we see that the Lie-algebra property of the 
Killing vectors may be expressed as 


Lele] = —£ele;] = Ía E ep. (2.59) 


Finally, in the phase space where by definition 0,p, = 0, we can also write the 
transformation (2.47) for the momentum as 
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OiPu = {J Pu = (Le; [pl),- (2.60) 


Thus we have established a close relationship between Poisson brackets and Lie 
derivatives, and between symmetries, Killing vectors and Lie algebras. 


2.5 Phase-space symmetries and conservation 
laws 


In the previous section it was established that continuous symmetries of man- 
ifolds are generated by Killing vectors and imply conservation laws for certain 
dynamical quantities (the generators J) along geodesics. However, not every con- 
stant of geodesic flow is necessarily connected with a Killing vector. This is most 
easily seen from the example of the hamiltonian itself, which is conserved but not 
of the form (2.45). Rather the hamiltonian is quadratic in the momenta. This 
leads us to consider the possibility of conservation laws for quantities which are 
higher order expressions in the momenta. 

Let J(z,p) be any function on the phase space of coordinates «“(\) and 
momenta p,(A), associated with a geodesic, which is non-singular in the momenta. 
In others words, we suppose it has an expansion of the type 


œ l 
J(x,p)= >> ii TY MME) Die Drug: (2.61) 
k=0 . 


Here the coefficients J can be taken to be completely symmetric in all k indices. 
We ask under what conditions such a function can be a constant of geodesic flow: 


dJ 
a = {H = 0. (2.62) 


To obtain the answer, insert the expansion (2.61) into this equation and compare 
terms at equal order of powers of momentum. At zeroth order one finds 


A, J =0, (2.63) 


i.e. J is a constant. Such a constant can always be absorbed in the definition 
of J and is of no consequence. From now on it will be dropped. 
At first order one finds back the Killing equation (2.44): 


DIM + DIM =0. 
These terms we have already discussed in detail. 


Finally, for k > 2 we find a generalization of the Killing equation for the 
co-efficients J) #--#* (a): 
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D,,, J” = 0, 2.64 
(u 


1 u2.Hk41) 

The parentheses denote full symmetrization over all component indices enclosed. 
The solutions of these generalised Killing equations are called Killing tensors. It 
is now obvious that the hamiltonian is included, because the metric gą» is co- 
variantly constant, and therefore automatically a solution of eq.(2.64) for k = 2. 
Moreover, this solution always exists, independent of the particular manifold con- 
sidered. In addition, for specific manifolds (specific metrics gay) other solutions 
may exist as well. For example, in spherically symmetric spaces the square of the 
angular momentum: 


= ; (2p — (a -p)}’) , (2.65) 


is conserved and quadratic in the momenta by construction. Other examples will 
be encountered later on. 

Just as the Killing vectors, the Killing tensors generate symmetry transforma- 
tions, not of the manifold but rather of the corresponding phase space spanned 
by the coordinates and the momenta (tangent vectors). In particular: 


ka 1 
da” = {J(æ, p) 2"} = -JO — > qg i TE O) PuPu (2.66) 
k=2 . 


We observe that the transformations generated by the higher-order Killing tensors 
are velocity dependent [18]. 

Finally we address the algebraic properties of the constants of geodesic flow. 
We have already remarked, that the Poisson brackets of two such constants yields 
another constant by the Jacobi identity (2.51). In general the rank of quantities 
on the left-hand side are additive minus one: 


LJO, JO} ~ JEH, (2.67) 


Now for | = 1 the rank of JEH- is the same as that of J). It follows, that 
(i) for the constants of geodesic flow linear in the momentum, corresponding to 
Killing vectors (k = 1 = 1), the algebra is closed in itself: it is a Lie-algebra G; 
(ii) for any fixed value of k, the brackets of a J™ and a J) close on the J™; 
however, the brackets of two elements J generates an element of order (2k—1). 
Hence in general the algebra (2.67) is not closed, unless one includes elements 
of arbitrary rank, i.e. all possible k. This implies that once higher-order invari- 
ants appear, one may expect an infinite series of them. Sometimes this happens 
indeed, as with the Virasoro algebra for the conformal transformation in two di- 
mensions [22]. However, in many cases the series stops because the higher order 
invariants are Casimir-type of operators, or they become trivial in that they are 
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just products of commuting lower-order elements (commuting here is meant in 
the sense of vanishing Poisson brackets). 

The Lie algebra G of Killing vectors is characterized by its structure constants 
Si; k. As the Poisson brackets of J) with all elements of the Lie algebra spanned 


by the J® are again of order k, one may choose a basis {elena for the 
Killing tensors of order k, labeled by the index A, and decompose the right-hand 
side of the Poisson bracket again in terms of the corresponding conserved charges 


JW. 


{I JP, = gP Ip). (2.68) 


This defines a set of higher-order structure constants. The Jacobi identities then 
imply, that the matrices (T;) P = —g;4” define a representation of the Lie algebra 
(2.53): 


GA Gio” = Jia Iie" = fy E Ipa” = [T Ti] 4 = fy; E Tpi. (2.69) 


Therefore the constants of geodesic flow 


JE (ap) = DO (0) pu Pugs (2.70) 


span a representation space of the Lie algebra G for the representation {T;}. 
The discussion of symmetries shows that group theory and Lie algebras can be 
important tools in the analysis of geodesic motion and the structure of manifolds. 


2.6 Example: the rigid rotor 


The above concepts and procedures can be illustrated by the simple example of 
a rigid rotor, which has a physical interpretation as a model for the low-energy 
behaviour of diatomic molecules. As such the example also serves to emphasize 
the usefulness of geometric methods in physics outside the context of general 
relativity. 

Consider two particles of mass mı and məz, interacting through a central po- 
tential V(r) depending only on the relative distance r = |r — rı|. The lagrangian 
for this system is 


L= Tia + a ~ V(r). (2.71) 
Defining as usual the total mass M = mı + my, and the reduced mass u = 


m m2/M, the lagrangian is separable in the center of mass coordinates defined 


by 


MR = myPry + MP2, T = T2 — Ti. (2.72) 
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In terms of (A,r) the expression (2.71) becomes 


L= it + ce ~ V(r). (2.73) 
As the center of mass moves like a free particle of mass M, it is convenient to 
work in the rest frame of the center of mass: R = 0. Then we are only left with 
the problem of describing the relative motion of the masses, described by the 
three coordinates represented by the vector r. 
Now suppose that the potential V(r) has a minimum for some separation 
r = ro, and rises steeply for all values of r near this minimum. Then the first 
excited vibration state of the molecule has an energy well above the ground 
state, and at low temperatures the only degrees of freedom that play a role in 
the dynamics are rotations of the molecule at fixed interatomic distance rọ. Thus 
the distance is frozen out as one of the dynamical degrees of freedom, and we are 
effectively left with only two coordinates: the angles (9, p) describing the relative 
orientation of the two masses; the potential V(r) = V(ro) is constant and may 
be ignored. This leaves as the effective action of the system 


Less = 5r (0? + sin? 0g). (2.74) 
This is the lagrangian of a particle of mass u moving on a spherical surface of 
radius ro. In mathematical physicists’ parlance the effective low energy degrees 
of freedom (6, Y) are sometimes called the modular parameters or moduli of the 
system, and the sphere representing the effective low-energy configuration space 
is called the moduli space. 

Up to an overall scale factor ur, the lagrangian Lepp represents an action of 
the form (2.35) in a 2-dimensional space with euclidean signature, described by 
the coordinates x! = 0, x? = y and metric 


1 0 
Jij = ( 0 sin2d (2.75) 


As the geodesics on a sphere are great circles, the great circles provide the general 
solution to the equations of motion. The spherical symmetry implies that the 
system is invariant under 3-dimensional rotations (rotations about each of three 
independent body axes of the sphere) and hence there are three Killing vectors 
generating these rotations. The corresponding constants of motion are the 3 
components of angular momentum L = (14, L2, L3), given explicitly by 


Lı = —siny pọ — cot b cos Y py, 
Tyg = cosg pọ — cot 0 sin Y py, (2.76) 
La = Po- 


29 


They span the Lie-algebra so(3) of the 3-dimensional rotation group: 


{Li, Lj} = Eijk Lee (2.77) 


As the Poisson brackets of L with the hamiltonian vanish, the hamiltonian can 
depend only on the Casimir constant of the rotation group, i.e. the total angular 
momentum squared L’. A direct computation shows, that up to the constant 
potential the effective hamiltonian is in fact proportional to the Casimir invariant: 


L? 
z: 


Hete = — 
H 215 


(2.78) 
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Chapter 3 


Dynamics of space-time 


3.1 Classical solutions of the gravitational field 
equations 


The Einstein equations (1.17) describe the gravitational fields generated by a 
specific distribution of material sources, as specified by the energy-momentum 
tensor T». In two or three dimensions, these fields carry at best a finite number 
of physical degrees of freedom, related to the sources and/or the topology [23]. In 
four or higher-dimensional space-time the gravitational field is a fully dynamical 
system with infinitely many degrees of freedom. In particular, dynamical gravi- 
tational fields can exist in the absence of material sources, in empty space-time. 
In the geometrical interpretation of general relativity such fields represent the 
intrinsic dynamics of space-time itself. 

In eq.(1.12) the curvature tensor was defined in terms of the parallel transport 
of a vector around a closed curve. A manifold is called flat if after parallel 
transport around an arbitrary closed curve in the manifold the image of any 
vector coincides with the original. Therefore a necessary and sufficient condition 
for a manifold to be flat near a given point is that the components of the Riemann 
tensor R „ò vanish there. This condition is more general than to require that 


LYK 
the metric takes the form 


ds? = +dz? + 5 dz?, (3.1) 


i=1 


with the sign depending on the euclidean or lorentzian signature of the manifold. 
Indeed, in any differentiable manifold one can always find a local coordinate 
system in the neighborhood of a given point, in which the line element can be 
diagonalized to this form, but this is not generally the case globally. Even in 
cases in which it can be done globally, eq.(3.1) is not the most general solution 
of the Einstein equations for flat empty space-time, because it is still allowed to 
perform an arbitrary local change of co-ordinates in this metric. The flatness- 
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criterion introduced above, amounting to the vanishing of the Riemann-curvature, 
is more useful as it holds in any co-ordinate system and in any topology. 

The existence of non-trivial solutions of the Einstein equations with vanishing 
cosmological constant in matter-free regions of space-time is possible because 
vanishing of the Ricci-tensor does not imply a vanishing of all components of the 
Riemann-tensor. A formal mathematical criterion for distinguishing solutions 
of the Einstein equations in empty space from those with material sources is 
through the Weyl curvature tensor, defined as the completely traceless part of 
the Riemann tensor: 


1 
Wuvraà = Raver = ) (Gun Ror — Jur Run = guà Pok + WrRux) 


(d? 
(3.2) 
+ TERTE (Jurga — JraJvr) R. 


Here d is the dimensionality of the manifold. If we similarly define W,» as the 
traceless part of the Ricci tensor: 


1 
Wav = Rav — goet Wav = Wins W,“ =0, (3.3) 
then we obtain a complete decomposition of the Riemann tensor in terms of 
traceless components and the Riemann scalar: 


Raver = Wvicd + Wi = Jur Waa = Gur Warr + WrW aur) 


1 

(d—2) (Sur 
(3.4) 
1 
+ daD (Jur Jvà — Imgur) R 

This decomposition allows a classification of the solutions of the Einstein equa- 
tions according to which of the independent components (Wava, Wuv, R) vanish 
and which don’t. In particular, metrics for which only the Weyl tensor is non- 
vanishing are solutions of the source-free Einstein equations Ray = 0. On the 
other hand, there also exist manifolds for which the only non-zero component 
is the Riemann scalar R, whilst the Weyl tensor and traceless Ricci tensor van- 
ish: the n-spheres r? + r? +... + x = p’, and all two-dimensional surfaces. 
For the n-spheres this follows essentially from the spherical symmetry, which im- 
plies the absence of any prefered direction, in combination with the permutation 
symmetries between the components of the Riemann curvature tensor. 

For two-dimensional manifolds, it is always possible to write the Riemann 
tensor as 
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1 1 
Raver = gI Ew env R = 5 (Gun Jud = GurJur) R. (3.5) 


For the (torsion-free) riemannian and pseudo-riemannian manifolds we consider 
in higher dimensions (d > 3), the number of independent components of the 
Riemann tensor R vrà, taking account of all its symmetries, is given by: 


N[Rwrrl = or) (3.6) 


This is counted as follows: the Riemann tensor is anti-symmetric in each of the 
index pairs [uv] and [KA]. Hence Rj,,j[<\) may be considered as a square matrix 
of dimension n = d(d —1)/1. Also as an (n x n)-matrix it is symmetric under 
interchange of these index pairs. As such it has n(n +1)/2 independent elements. 
Finally the cyclic symmetry in the first three (lower) indices of the fully covariant 
components of the Riemann tensor expressed by the Bianchi identity (1.15) adds 
d(d — 1)(d — 2)(d — 3)/4! algebraic constraints, eliminating an equal number of 


independent components. Combining these results one finds precisely the number 
N[Ruvra] in eq.(3.6) above. 

For d > 3 the number of components of the traceless Ricci tensor and Riemann 
scalar is 

d(d+1 
NW] = ach) — 1. N[R] = 1. (3.7) 

(for d = 2 the Riemann curvature tensor has only one independent component, 
corresponding to the Riemann scalar). Therefore for d > 3 the number of inde- 
pendent components of the Weyl tensor is 


N[Wyma] = Zda -1)(d+2)(d—3). (3.8) 


Clearly this vanishes for d = 3, and is positive non-zero only for d > 4. Thus 


we find indeed that non-trivial dynamical solutions to the source-free Einstein 
equations can exist only in space-times of dimensions d > 4. As the minimal 
case d = 4 happens to be the number of dimensions of our (macroscopic) world, 
we observe that for this case the Weyl tensor and the Ricci tensor each have 10 
algebraically independent components. Therefore half of the components of the 
curvature are determined by the material sources, whilst the other half describes 
the intrinsic dynamics of space-time. 

In the geometrical interpretation of the gravitational field and its dynamics, 
Einstein’s famous principle of the equivalence of accelerations and local gravita- 
tional fields can be formulated as the hypothesis that space-time is a manifold, in 
each point of which one can establish a locally flat (cartesian) coordinate system 
of the form (3.1). Of course a fully geometric interpretation of all of physics, 
and also the equivalence principle, breaks down in those points or regions where 
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material sources are located!. But also solutions of the source-free Einstein equa- 
tions can violate the equivalence principle locally if there are singularities — for 
example points where curvature invariants become infinite and/or the number of 
dimensions changes. Physics at these points is no longer described, at least not 
completely, by the Einstein equations, and space-time is not necessarily a mani- 
fold in the strict sense. In some cases it can still be a manifold with boundaries, 
where specific boundary conditions hold; but these boundary conditions then 
represent external data and the gravitational field is not determined by general 
relativity alone. 

Nevertheless singular solutions are important because they are relevant for 
physics as we observe it in the universe: cosmological solutions of the Friedmann- 
Robinson-Walker type with a ‘Big Bang’-singularity, and black-hole solutions 
describing stellar collapse. The singularities in these cases are not necessarily 
physical, as the short-distance physics may be regulated by Planck-scale mod- 
ifications of the theory, which are presently unobservable. However that may 
be, here we simply accept the possibility of singularities and discuss examples 
of both non-singular and singular geometries arising as solutions of the Einstein 
equations. They illustrate many interesting aspects of gravitational physics. 


3.2 Plane fronted waves 


To show that in four dimensions non-trivial gravitational fields can indeed prop- 
agate through empty space-time, we now consider a class of solutions known in 
the literature as plane-fronted gravitational waves [4, 5]. This name derives from 
the property that there is a co-ordinate system in which these field configurations 
(i) move with the speed of light along a straight line; (72) have a flat planar wave 
front, and (iii) are of finite duration. When discussing the physical properties of 
these solutions it will become clear that these wave-like solutions actually behave 
more like dynamical domain walls. 

We construct a planar wave solution with the direction of motion along the 
x-axis, and the transverse directions defined by the (y, z)-plane. Thus we look 
for solutions of the Einstein equations of the form 


—1 


Jur(X,t) = Pet) . (3.9) 


The corresponding space-time line-element is 


l Attempts to explain matter as a manifestation of physics in higher-dimensional space- 
time, as in Kaluza-Klein models or the more sophisticated string theories, aim to redress this 
situation. 
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ds? = —c’dt? + dx? + f?(x,t)dy’? + g?(a,t)dz’. (3.10) 

We observe the following properties of these metrics: 
— in the (x, t)-plane space-time is flat (Minkowskian); 
— in the y- and z-directions space-time is generally not flat, but the value of the 
metric depends only on the distance along the z-axis and on time, not on the val- 
ues of y or z themselves. Hence in the co-ordinate system chosen the components 
of these gravitational potentials are the same everywhere in the (y, z)-plane. In 
particular they do not fall off to unity at spatial infinity, showing that the space- 
time is not asymptotically Minkowskian at large space-like distances. Of course, 
they share this property with plane electromagnetic waves, which at a fixed time 
also have the same amplitude everywhere in the transverse directions; 
— the metric has no inverse if f = 0 or g = 0; then the connections are singular. 
We must take special care when this happens. 

With the Ansatz (3.9) for the metric one can compute the components of the 
connection using eq.(1.9). Using the notation æ“ = (ct, x,y,z) for u = (0,1,2,3), 
the non-vanishing components have been assembled in the following list: 


(3.11) 


The next step in solving the Einstein equations is to compute the components of 
the Riemann curvature tensor using eq.(1.14). The only independent covariant 
components Raris = Rana Ion not identically zero are the following ones: 


1 
Ro202 = 2 F ft, 


Rioi2 = Í few 
i (3.12) 
Roz12 = P Í fts, 


ji 
R2323 = fa( Sege = z Fige). 
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Before discussing the general solution of the Einstein equations for this class of 
metrics, we first use the above results to answer the question which of the metrics 
singled out by the Ansatz (3.9) describe flat Minkowski space. If we require all 
components of the curvature tensor (3.12) to vanish independently, the most 
general solution for f(x,t) and g(a, t) is 


f(a,t) = a+ px + ret, 
glz, t) = b+ qx + let, (3.13) 
with the additional constraint 
pq =rl. (3.14) 


Thus metrics of the form (3.9), and corresponding line-elements (3.10), describe 
a flat space-time if f and g are linear functions of z and t, with the additional 
constraint (3.14) between their co-efficients. 


The components of the Ricci tensor R, now follow by taking a covariant trace 
as in eqs.(1.16). The non-vanishing ones are 


(3.15) 


In the absence of material sources Einstein’s equations require the Ricci tensor 
to vanish. Assuming f and g themselves not to vanish identically, some manipu- 
lation of these equations then leads to the following five necessary and sufficient 
conditions: 
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Jez Z fa = 0, 
Jes — 5 91 = 0, 
Se Ja — 5 fan = 0, (3-16) 
Ifo + fue = 9, 
Jfis + foe = 0. 


The first two equations are free wave equations in two-dimensional (x,t) space. 
They imply that f and g can be decomposed into left- and right-moving waves: 


flæt) = fyla tet) + f(x — c), 


(3.17) 


glz,t) = g4(z+ct) + g-(x -— ct), 
where (f_(u),g-(u)) and (f4 (v), g+}(v) are functions only of the single lightcone 
variable u = z — ct, and v = z + ct, respectively. Each of these components 
then satisfies the wave equation by itself, and at this stage they can be taken to 
be completely independent. The only ambiguity here is that a possible constant 
term can be divided in an arbitrary way between the left- and right-moving 
components of the solution. 

If we now substitute this decomposition of f and g in the third equation 
(3.16), and assume that the first derivatives of the components with respect to 
u or v (i.e. fi and gi) are all non-zero, then the only solutions of the Einstein 
equations correspond to flat Minkowski space. This can be seen by dividing the 
equation by gg- to obtain 


B E 


gh 7 gL 


where A on the right-hand side of the last equality is a finite, non-zero constant. 


=, (3.18) 


This last equality holds because the two previous expressions are functions of 
different variables, hence they can be identically equal only if they are constant. 
It then follows, that 


fi = Ags, fl = —dg!. (3.19) 
Inserting these results into the last two equations (3.16), addition and subtraction 
gives the equivalent conditions 


(àg + f) gf = 9, (àg— f) g” = 0. (3.20) 
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The solution gi = g” = 0 implies ff’ = f” = 0. This leads to linear dependence 
of f and g on x and t: they are of the form (3.13), (3.14) and correspond to flat 
Minkowski space. 

On the other hand, in case one of the second derivatives g¥ is non-zero, one of 
the co-efficients (Ag f) must vanish; note that they cannot vanish simultaneously 


for non-zero f and g. Hence eq.(3.20) gives either 


f = Ag and gy = ff =9, (3.21) 


or 


f = —Ag and g! = f" =0. (3.22) 


But combining this with the results (3.19), we find again that all second deriva- 
tives must vanish and that the geometry is that of flat space-time. 

We conclude that in order to get non-trivial solutions, at least one of the first 
derivatives fi or gi, must vanish. The last three equations (3.16) then imply 
that, after absorbing any constants terms in the remaining components, the only 
solutions are 


f(x,t) = fale + ct), glz, t) = g4(x + et), (3.23) 


or 


f(x,t) = f-(x — et), g(z,t) = g-(x — ct), (3.24) 


with the additional condition 


f” g" 

F + 7 0. (3.25) 
Although f and g represent waves traveling either to the right or to the left at the 
velocity of light, the non-linear character of the Einstein equations exemplified by 
the additional condition (3.25) do not allow us to make arbitrary superpositions 
of these solutions. Interestingly, viewed as a two-dimensional theory in (x, t)- 
space f and g represent chiral bosons. Wherever necessary in the following we 
chose to work with right-moving fields f(a,t) = f-(u) and g(#,t) = g_(u) for 
definiteness. 

An example of a simple solution of eq.(3.25) is [9]: 


f(a,t) = cosk(a— ct) = cos ku, 


g(x,t) = coshk(a — ct) = cosh ku. (3.26) 


Note that at u = (2n+1)m/2k the component f has a zero, whilst the solution for 
g(x, t) grows indefinitely with time or distance. However, this does neither imply 
that the geometry is singular, nor that the gravitational fields become arbitrarily 
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strong: the apparent problems turn out to be coordinate artifacts. This point is 
discussed in detail in the next section. 

Let us consider a wave which reaches x = 0 at time t = 0, and extends to 
x = —L, i.e. the solution (3.26) holds for —L < u < 0. If |kL| < 2/2 the metric 
components f? and g? are strictly positive. For u > 0 space-time is flat; we can 
choose a co-ordinate system in which f = g = 1. Note that the wave solution 
(3.26) assumes these same values at u = 0, whilst the first derivatives match as 
well. As a result the solutions are smoothly connected and satisfy the Einstein 
equations even at the boundary u = 0. 

For u < —L space-time is again flat, but now we have to use the general 
solution (3.13): 


f(z,1) = a+ p(z- ct), 
b + q(x — ct), (x —ct) < —L. (3.27) 


Q 

& 
~ 
II 


The co-efficients (a,b; p,q) in this domain of (x,t) values are to be determined 
again by requiring the metric to be continuous and differentiable at the boundary 
u = —L between the different domains of the solution. This condition is sufficient 
to guarantee that the Einstein equations are well-defined and satisfied everywhere. 
Imposing these requirements we find 


coskL = a-— pL, ksinkL = p 
(3.28) 
coshkL = 6b —qL, —ksinhkL = q. 
For f and g in the domain (æ — ct) < —L we obtain the result 
f(z,t) = coskL + k(L+e«-—ct)sinkL, 
(3.29) 


glx,t) = coshkL — k(L +z- ct) sinh kL. 


Conclusion: we have constructed a solution of the Einstein equations in empty 
space which interpolates between two Minkowski half-spaces, connected by a fi- 
nite region of non-zero curvature, where the gravitational fields have non-zero 
gradients and test particles are accelerated with respect to the initial inertial 
frame. After passage of the wave the test particles have acquired a finite, but 
constant velocity. This shows that they have been accelerated, but the accelera- 
tion has ceased after t > (x + L)/c. 


The plane-fronted wave constructed here is a well-defined solution of the equation 
Ra = 0 everywhere, including the wave fronts x = ct and z = ct — L, as the 
contributions of the metric components g,, = f? and g., = g° always compensate 
each other in eq.(3.25). However the Riemann curvature itself, which coincides 
with the Weyl curvature in this case, does not vanish everywhere. Indeed, we 


find from eq.(3.12) 
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A 
Ro202 = Riz12 = + Ro212 = ff ’ 


Rosos = Rizizs = Rosis = gg”, (3.30) 


R2323 = 0. 


Inserting the explicit solution for f and g for each of the three regions x > ct, 
—-L<a«-—ct<Oand x < ct— L one finds that the curvature is zero in the two 
Minkowski half-spaces before and after the wave, as expected, but finite and non- 
zero in the interior of the wave. Indeed, the second derivatives f” and g” have a 
finite discontinuity on the wave front of magnitude =k? at x = ct, and —k? cos k L, 
resp. k? cosh kL, at x = ct — L. This finite discontinuity then also appears in the 


curvature components, where it must be attributed to a finite discontinuity in 
the Weyl curvature, and thus can not be attributed to a (singular) distribution of 
sources. But the Riemann-Christoffel connections, which depend only on the first 
derivatives of the metric, are continuous. Then the geodesics and the motions 
of test particles in the space-time described by the wave are continuous, and it 
follows that the jump in the Weyl curvature presents no physical inconsistencies. 


3.3 Nature of the space-time 


The planar wave constructed above can be characterised as a region of non- 
vanishing curvature sandwiched between two flat space-times and propagating 
with the speed of light —either in the positive or the negative x-direction— such 
that one of the two flat space-time regions grows at the expense of the other. 
In this section I establish the relation between the two flat space-times, before 
and after the wave. The analysis shows that they are not related by a simple 
Lorentz transformation: a set of equidistant test particles carrying synchronised 
ideal clocks and thus defining an inertial co-ordinate system before the wave 
are neither equidistant nor synchronised in any inertial co-ordinate system after 
the wave. Therefore the two Minkowski space-times on either side of the region 
of non-zero curvature are inequivalent, in the sense that the Lorentz group of 
rigid transformations leaving the inproduct of four-vectors invariant in these flat 
regions are embedded differently in the group of local co-ordinate transformations 
which constitutes the symmetry of the full theory. 

This situation is also encountered in scalar field theories with broken sym- 
metry, when there is a boundary separating two domains with different vacuum 
expectation values of the field. One can even show that in some models these 
two forms of spontaneous symmetry breaking —of a scalar field and of the grav- 
itational field— are connected [24]. Therefore we find a new interpretation for 
the planar gravitational waves: they act as domain boundaries separating differ- 
ent, inequivalent flat-space solutions (‘vacua’) of the Einstein equations. These 
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domain boundaries are dynamical: they are not static but move at the speed of 
light; for this reason they are sometimes refered to as ‘shock waves’. However, 
the planar wave solutions discussed here have finite extent, whereas usually the 
term shock-wave is reserved for waves with a delta-function profile. 

Regions of space-time enclosed by such dynamical boundaries clearly grow 
or shrink at the speed of light. If we consider our planar wave solution as an 
approximation to a more general domain-wall geometry with a radius of curvature 
large compared to the characteristic scale of the region of space-time we are 
interested in, then we can use it to describe the dynamics of gravitational domain 
structure, which may have interesting implications for example for the physics of 
the early universe. 

Having sketched this general picture, I now turn to describe the sandwich 
structure of the planar wave geometry more precisely. A convenient procedure 
is to introduce new co-ordinates X” = (cT, X,Y, Z) related to the original z“ = 


(ct, x,y,z) by 


A 
cT = cd—-, 
2 
A 
X = zr-7, 
2 (3.31) 
Y = fy, 
Z = gz, 
where in terms of the original coordinates x 
Aa") = y ff + qq. (3.32) 


This coordinate transformation takes an even simpler form when expressed in 
light-come variables 


U = X-cP =a-cd =u, 
(3.33) 
V = X+ =eric-A=v-A. 
Indeed the line-element (3.10) now becomes 
d? = —e?dt? + dx? + Pdy + Pd? 
(3.34) 
= dUdV + K(U,Y, Z)dU? + dY? + d2?, 
with 
H g" g" 
K(U,Y, Z) == Y° 452 => (Z-Y"). (3.35) 


f g g 


4? 


Both for t < x/c and t > (x + L)/c we have f” = g” = 0, and therefore the line 
element (3.34) is manifestly flat Minkowskian on both sides of the planar wave. 
However, inside our particular wave solution one has 


H i 
Lo oop (3.36) 


f g 

Then K = k?(Z? — Y°), independent of U. The coordinate system X” with 
the line-element (3.34) is the closest one to an inertial frame at all times for an 
observer at rest in the origin one can construct. Indeed, the point x = y= z =0 
corresponds to X = Y = Z = 0 at all times, and both before and after the 
wave the X” define a global inertial frame. Furthermore, inside the wave the 
new coordinate system deviates from an inertial frame only for large values of 
(Z? —Y?). We also observe that in this coordinate system the line element (3.34) 
never becomes singular, showing the previous zero’s of det g, corresponding to 
zero’s of f(u), to be coordinate artifacts indeed. 

Computing the various geometrical quantities: metric, connections and cur- 
vature components, in the new coordinates produces the following summary of 
results: 

— The only non-vanishing components of the inverse metric are: 


gY = gh! =2, gY =4k(Y? — Z’), 


(3.37) 
gY = 97% =1. 
— The only non-vanishing components of the connection are: 
Tyg =RY, Tye = Poy = —2k°Y, 
(3.38) 
— The only non-vanishing covariant components of the curvature are 
Ruyvy = —Ruguz = —k’. (3.39) 


Thus in these coordinates the curvature components are constant inside the wave, 
and we conclude that our particular planar wave solution is actually like a block- 
wave: at the boundaries U = 0, U = —L the curvature components jump from 
zero to the values (3.39) and back to zero, whilst in between they are constant. 
This makes it easy to visualize how our particular solution (3.26) can be used to 
construct general planar wave solutions with polarization along the Y- and Z- 
axes: essentially it amounts to glueing together sufficiently narrow block waves 
of fixed curvature to obtain solutions in which the curvature components vary 
according to some arbitrary profile. 

The result (3.39) also allows one to verify once more by a one-line computation 
that the Ricci tensor vanishes identically for these planar-wave geometries. 
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With these results in hand, our next aim is to establish the relation between the 
flat half-spaces before and after the wave. We have already observed that the 
origin « = y = z = 0 becomes the origin X = Y = Z = 0 for all times ¢ or 
T, respectively. Next observe, that the coordinate transformation (3.31) is not 
a Lorentz transformation: the coordinates (t, x,y,z), defining an inertial frame 
before the passage of the wave, represent a non-inertial frame after passage. As 
a check, note that the jacobian of the transformation is 


OX" 


Ox” 


Indeed, a system of synchronized, equidistant clocks at rest in the first coordinate 


= fg Fl, foru = z — ct £0. (3.40) 


system, which operationally defines an inertial frame before the wave arrives, is 
neither synchronized nor equidistant in the final inertial frame. Moreover, in 
general it cannot be made so by a simple Lorentz transformation: the clocks have 
been accelerated with respect to each other by the wave, and ultimately define 
a non-inertial frame at late times. This is verified by the explicit calculation of 
geodesics in the next section. 

¿From this analysis it will be clear, that in principle it is possible to deter- 
mine the passage, and measure the characteristics, of a planar gravitational wave 
by measuring the acceleration it imparts to test masses in the laboratory. In- 
deed, such measurements can determine the functions (A, f, g) from the relation 
between (t, x,y,z) and (7, X,Y,Z). Although in this form the statement holds 
only for waves of the type discussed, the principles for other types of waves are 
the same, even if the functions characterizing the relation between the initial and 
final Minkowski spaces are generally different. 


3.4 Scattering of test particles 


In order to study the response of test masses to passing gravitational waves, 
it is necessary to solve for the geodesics of space-time in the gravitational field 
generated by the wave. We are particularly interested in the worldlines of particles 
initially at rest in a specific inertial frame (the laboratory). Observing frame gives 
the scattering data, from which the properties of the gravitational wave can then 
be inferred relatively easily. 

In principle it is very easy to find geodesics: they correspond to the world- 
lines of particles in free fall. Consider therefore a particle at rest in the initial 
Minkowski frame; its four co-ordinates are 


z”(T) = (cr, £1, &, &), (3.41) 


where the three vector € = (£1, £2, 3) fixing the position of the particle is constant; 
note also that by definition for a particle at rest in an inertial frame the co- 
ordinate time equals the proper time T. 
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In the co-ordinate system we use, a particle with co-ordinates (3.41) is actually 
always in free fall, its co-ordinates in the x“-frame do not change if it enters a 
non-trivial gravitational field. Of course, in an arbitrary gravitational field the 
x"-frame is not an inertial frame, and even if at late times space-time is flat 
again the particle will have acquired a finite, generally non-vanishing four velocity. 
However, its position in the final Minkowski frame can be found straightforwardly 
by applying the co-ordinate transformation (3.31). 

Let us first confirm the above statements by a simple mathematical argument. 
Then we turn to apply the result to find the actual motion of the test particle 
in the final Minkowski frame. We begin by noting that in the z“-frame the 
metric for the gravitational wave is always of the form (3.9), and therefore the 
only non-vanishing components of the connection are those given in eq.(3.11). In 
particular, in the z“-frame the components loo“ always vanish. 

On the other hand, for a particle with the worldline (3.41), the four velocity 
is purely time-like and constant: 

dx" 


u“ = —— = (¢,0,0,0). (3.42) 
T 


Then the four acceleration in this frame vanishes identically, as does the gravita- 
tional force: 


Dy“ uu’ = 0. (3.43) 


As a result the geodesic equation (1.10) is satisfied, proving that (3.41) represents 
the worldline of a particle at all times. 


What we are particularly interested in, is the position and velocity of the parti- 
cle at the moment it comes out of the gravitational wave. It has already been 
remarked that the origin € = 0 remains at rest and coincides with the origin 
of the final Minkowski frame. Therefore at late times the origin has laboratory 
coordinates 


Xë (r) = (er, 0,0,0). (3.44) 


Let us compare this with the four co-ordinates of a particle away from the origin 
in the (y, z)-plane: € = (0, £2, €3) # 0. For such a particle the values of the metric 
components at late times are 


f(O,7) = a — pet, 
nn — b «cr. cT > L. (3.45) 


Here (a,b; p,q) have the values quoted before, eq.(3.28). For A one finds in the 
same domain ct > L: 


A = ap + ba — (PÉ + PG) cr. (3.46) 
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As a result, the Minkowski co-odinates of this particle at late proper times are 


Xi (7) = (cTi, Xi, Yi, Za), (3.47) 
with 
1 CT 
chi = -3 (apé + b48) + D (24+ P6476), 
o | 2 2 CT ò 949 22 

X = -7 (apé T baé?) T > (p & +q E), (3.48) 
Yi = aks = p&ocT, 
Z, = bs — g&scr. 


Note, that clocks which were originally synchronized —same value of t = 7 at 
early proper time— are no longer so after passage of the gravitational wave: 7} is 
not equal to To = 7 at late proper time. Before we can calculate the position of 
the test particle at synchronized co-ordinate time T = 7, in the final Minkowski 
frame, we first have to correct for this gravitationally induced position-dependent 
time shift. 

To determine the instantaneous proper distance between the particle with 
worldine Xj{‘(7) and the origin, we have to ask at what proper time 7 the time 
co-ordinate T; equals the time Tọ measured by a clock at rest in the origin: 


Ti(7) = 7, (3.49) 
with the result 


2 2 2 
z = Zor t apt, + bats (3.50) 
2+ pE +E 


The position at this proper time, X“ = Xj'(7), can now be evaluated: 
7 = (apts + 6963) + cT (pE + g6) 


X, = ee ae a 
' 24+ p?& + 9&3 


7 (aq — bp) g&3 +2 (a — pen) 
n = ra &1?|? 3.51 
15N ( 2+ PG +P (3.51) 
> — (aq — bp) pf + 2(b— ¢er) 
Z =% a 
2 +p E +E 
Finally we can also determine the velocity of the test particle: 
v 1 242 | 2¢2 
e 24 PEHEE (p E3 + 4'63, —2pha, 2463) , (3.52) 
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and as a result the kinetic energy is 
ly = me PEED UHE +P) 
2 (2+ PE + P&) 


From this analysis we draw the conclusion, that not only in general test particles 


(3.53) 


are accelerated by the gravitational wave, but also that the acceleration and 
hence the final velocity depend on the position of the particle. Therefore no 
Lorentz transformation with any fixed boost velocity v can bring all test particles 
initially at rest w.r.t. one another back to rest simultaneously after passage of 
the wave. This result explicitly demonstrates the inequivalence of the initial and 
final Minkowski space-times. 


3.5 Symmetry breaking as a source of gravita- 
tional waves 


The planar gravitational wave solutions described in this chapter arise as transi- 
tion regions between flat domains of space-time, which are related by a coordi- 
nate transformation that cannot be reduced to a (linear) Lorentz transformation. 
We have therefore interpreted them as domain boundaries between regions of 
space-time trapped in inequivalent classical vacua of the gravitational field (in 
the absence of a cosmological constant). This is possible because flat space (zero 
curvature) is associated with a value of the gravitational field, the metric guv, 
which in an appropriate coordinate system (gauge) is constant, non-vanishing 
and Lorentz-invariant, but not invariant under general coordinate transforma- 
tions. In this sense the general coordinate invariance of the Einstein equations 
is broken spontaneously by the vacuum solutions corresponding to Minkowski 
space-time. 

A domain structure arising from the presence of degenerate but inequivalent 
classical vacua is a well-known consequence of spontaneous symmetry breaking in 
scalar field theories. In this section this point is elaborated on and it is shown, that 
dynamical domain boundaries in scalar field theories (scalar-type shock waves, 
if one prefers) can actually be sources of planar gravitational waves of the kind 
analysed in detail above [24]. 

The simplest scalar field theory with a continuously degenerate set of classical 
vacua is that of a complex scalar field with a ‘mexican hat’ type of potential, as 
described by the lagrangian 


a 
vg 


In this lagrangian cà is a dimensionless coupling constant and u is an inverse 


À 2\? 
= 99,006 3 (wo), (3.54) 


length (wh/c is a mass). In the natural system of units in which h = c = 1 we can 
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take A to be dimensionless and u to represent a mass. In the following equations 
one then can substitute c = 1 everywhere. 

This lagrangian (3.54) is invariant under rigid U(1) phase shifts ¢ > ġe”, 
with 7 constant. As a result for u? > 0 there is an infinitely degenerate set 
of solutions of the equations of motion minimizing the potential (refered to as 
classical vacua), of the form 


d= te”, (3.55) 


with @ an arbitrary real constant. Clearly the solutions themselves are not in- 
variant under the U(1) transformations: the classical vacuum solutions break the 
U(1) symmetry spontaneously. However, because of the U(1) symmetry of the 
dynamics these solutions are physically indistinguishable, and the absolute value 
of # is unobservable. On the other hand, variations in the vacuum angle between 
different regions of space-time can have observable effects [26, 27]. 

Parametrizing the complex scalar field as 6(x) = p(x) /\/2, the lagrangian 
becomes 


1 1 À 2P \ 1 
Jaa’ = ZI Oupðvp _ 8 ( 2L ac) — zP 9“ 0,080. (3.56) 


This shows that in the spontaneously broken mode of the U(1) symmetry the 
theory desribes actually two real scalar fields: the field y = p — u,/2/À, obtained 


by shifting p by its classical vacuum value u/2/A, with a mass (or inverse co- 
herence length) equal to u, which can become arbitrarily large; and the strictly 
massless field o = noJ2/, representing the Goldstone excitations arising from 
spontaneous symmetry breaking [26, 27]. 
Clearly at energies small compared to the scale u set by the mass of the field 
y the excitations of this massive field play no role and one only has to consider 
the Goldstone field o, described by the effective lagrangian 
=l e= L ry Lgro,oð 3.57 
Vag rG T gI o (3.57) 
where we have added the Einstein-Hilbert action for the gravitational field as 
this also describes a low-energy degree of freedom. Note that with the potential 
defined in eqs.(3.54), (3.56) there is no cosmological constant in the presence of 
the vacuum value of the scalar field, p = u/2/A or equivalently y = 0. 
The field equations derived from this lagrangian are the Klein-Gordon equa- 
tion for the Goldstone field oø: 


1 
g = — d, y —gg 0o = 0. 3.58 
=a, (3.58) 
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and the Einstein-equations with a source term from the energy-momentum of the 


Goldstone field: 


8SrG 
Ryu = a ð, ,o0,0. (3.59) 
Remarkably this system of equations admits exact planar wave solutions similar 
to the free (sourceless) Einstein equations. Using the same Ansatz (3.10) for 
the metric we find running planar gravitational waves (f(a + ct), g+(x + ct)) 


in combination with a scalar wave ox(a + ct) running in the same direction. 


The continuity conditions require that this scalar wave interpolates continuously 
between a domain in which o = go before the wave and a domain where o = g1 
after the wave, where both (oo, 01) are constants. 

The simplest solution of this kind is a profile of o interpolating linearly be- 
tween these values: 


o_(u) = oo, for u > 0; 
o_(u) = oo — kupey/2/d, for — L <u <0; (3.60) 


o_(u) = oi = oot kLyy/2/X, foru < —L, 


with similar expressions for the left-moving solutions o,(v). It follows that the 
slope parameter k in (3.60) is given by 


01 — 09 À 
k = =. 3.61 
z \? (3.61) 
At this point it is useful to recall that o actually represents the phase of the 
complex scalar field ¢. Hence with a fixed value of the modulus p a monotonic 
linear increase of o implies a pure monochromatic oscillation of the real and 


imaginary components of the scalar field: 
_ iku _ !  _ ik(w—ct) 
plu) = doe" = -= e (3.62) 


where the last equality assumes p to equal its classical vacuum value and the 
gauge choice og = 0. It follows that the slope parameter k actually represents 
the wavenumber of these ‘large’ oscillations of @ interpolating between different 
classical vacuum values. 

It remains to solve the Einstein equations (3.59). For the right-moving wave 
solution (a_(u), f_(w), g_-(w)) these equations reduce to a slightly modified ver- 
sion of eqs.(3.16), the only modification being a constant source term in eq.(3.25): 


— + = hk’, for -L<u<0. (3.63) 


f" g" l67rG? 
f g ctr 
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Outside the region —L < u < 0 the right hand side vanishes, and we have the 
trivial flat-space solutions f” = g” = 0. 

It is not very difficult to establish the existence of solutions to equation (3.63). 
In fact there are more solutions than the purely gravitational case, because first 
there are the planar gravitational block waves of the type 


f(u) = cos kiu, glu) = cosh kau, (3.64) 


with the constraint between the parameters (%1, K2): 


l6nGy? 

2 2 2 

ki — ky = —— k. 3.65 
1 2 ct) ( ) 
These waves have a purely gravitational contribution from a non-vanishing Weyl 
tensor. But in addition there are also purely oscillatory solutions 


f(u) = cos kyu, glu) = cos kou, (3.66) 


with a relation between the wavenumbers (k1, K2) of the type 


2 


It is quite clear that these purely oscillatory waves require a non-vanishing of 
the Ricci tensor: they can exist only in the presence of the scalar field o_(w), 
because in the limit k — 0 one finds simultaneously that both xı, — 0. This is 
in contrast to the hyperbolic type of solution (3.65), where k = 0 only implies 
Kk, = Kg. Therefore one might think of the first type of solution as radiative 
waves, which have kept the quadrupole character known from the propagation of 
small perturbations in the linearized version of gravitational field theory, whereas 
the second type are matter-induced waves with a dipole-like behaviour. 

By performing the coordinate transformation (3.31) we can again construct a 
reference frame which is manifestly Minkowskian before and after the wave, and 
has constant curvature components inside the wave. This is based on the choice of 
the simple solution (3.60) for the scalar field, with constant slope parameter k, i.e. 
a monochromatic @-wave. As in the case of the purely gravitational waves, more 
general solutions with smooth behaviour at the wavefront can be contructed. 

The main lesson to be learnt from the results in this section is that the in- 
terpretation of the gravitational waves as moving domain boundaries between 
inequivalent classical vacua has its parallel in the case of scalar fields with bro- 
ken symmetry, and that these scalar waves are actually accompanied by planar 
gravitational waves of the type constructed above. 
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3.6 Coupling to the electro-magnetic field 


The coupling of planar gravitational waves to scalar fields of the Goldstone type 
is the result of the existence of field configurations with constant slope (in do- 
main boundaries), for which the energy-momentum tensor has constant non-zero 
components in half of the directions of space-time. A similar type of energy- 
momentum density can be generated by constant electric and magnetic fields 
and again one finds that these are accompanied by gravitational waves [25]. 

We consider the simplest case of constant electric and magnetic field perpen- 
dicular to each other, for example in the (y, z)-plane: 


> > 


E = (0, E,0), B = (0,0, B), (3.68) 


and a metric which is non-flat only in the (zz) direction: 


d? = —c'dt? + dx? + dy? + gdz’, (3.69) 
These fields provide a solution of the combined Einstein and Maxwell equations: 


8rG 8TEoG 1 
Rmuv = -I Tuv = — 2 (Fat, — IE?) 5 (3.70) 


and 


D,F“” = 0, (3.71) 


provided the electric and magnetic field are constant and equal modulo c in 
magnitude: 


E = +cB = constant, (3.72) 


whilst the metric component g., represents a planar wave 


1 327 E0G 
g(x + ct) = cos ze + ct), k = “ 


F°. (3.73) 


C 


Here the sign in the argument (i.e. left or right moving wave) depends on the 
sign in (3.72). 

It is straightforward to check that this is a solution of the Einstein equation 
by substituting the g as above and f = 1 in eq.(3.15) and writing out the energy- 
momentum tensor of the electro-magnetic field: 
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Tw = 
I(E? + cB?) cBE 0 0 
cBE HE? + eB?) 0 0 
Ege 
0 0 —3(E? — cB?) 0 
0 0 0 “(B? — êB?) 
(3.74) 
Observe, that the metric component g, itself becomes 
ə 1l 
G2 =P =3 (1+cosk(a+ct)), (3.75) 


ans oscillates with frequency 


3E 
w = ke = — VTEC. (3.76) 


It may seem surprising and worrying that constant electric and magnetic field 
generate gravitational waves of such large amplitudes; however, note that in 
laboratory units the frequency is 


E 
[IV /m] 


For ordinary electric fields in the laboratory this means that it takes a full lifetime 
of the universe or longer to go through one complete cycle. Even at the high- 


v=1.3 x 107" Hz. (3.77) 


est fields imaginable in theory, when the vacuum breaks down due to electron- 
positron pair creation: 


TMC? 
CA, 


where e is the electron charge, me its mass and Àe its Compton wave length, the 
frequency can not get higher than about 1 Hz. Therefore the effect will not be 
measurable under ordinary circumstances. 

Nevertheless it does seem remarkable that constant electro-magnetic fields 
generate time-dependent (oscillating) gravitational fields. This effect, as well 


Berit = = 4.1 x 10! V/m, (3.78) 


as the scalar-generated waves discussed in the previous section, has quite close 
similarities to the Josephson effect in superconductors, where a constant potential 
generates an oscillating current. In that case the current is related to the change 
in the phase of the electron condensate, which is mathematically analogous to 
the phase of the complex scalar field represented by the Goldstone mode o. 
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Chapter 4 


Black holes 


4.1 Horizons 


The example of gravitational shock waves/domain walls has shown that space- 
time in four dimensions is a truely dynamical arena for physics even in the absence 
of matter, although adding matter —for example in the form of scalar or elec- 
tromagnetic fields— makes the possible structures even richer. In this chapter 
we turn to another class of very interesting structures in four-dimensional space- 
time: static or quasi-static solutions of the Einstein equations which represent 
gravitating extended bodies, whose static fields become strong enough to cap- 
ture permanently anything that gets sufficiently close to the central core. As 
these quasi-static objects even capture light they are called black holes, but it 
has turned out that their properties are much more interesting and peculiar than 
this ominous but dull sounding name suggests. 

Black holes possess two characteristic features. The first, which has given 
them their name, is the existence of a horizon: a closed surface through which 
objects can fall without ever being able to return. The second is the presence of 
a singularity, a locus of points inside the horizon where the curvature becomes 
infinite. It has even been conjectured by some authors that these two properties 
are so closely linked that any space-time singularity should always be hidden 
behind a horizon. This conjecture bears the name of the Cosmic Censorship 
hypothesis. Although as an unqualified mathematical theorem the hypothesis 
is certainly not correct, it could be true for singularities that arise in realistic 
processes from the gravitational collapse of macroscopic bodies, like imploding 
massive stars. 

In this chapter we restrict our attention to the special black-hole solutions that 
might be called classical ‘eternal’ black holes, in the sense that at least outside 
the horizon they describe stationary gravitational fields, and do not require for 
the description of their behaviour the inclusion of their formation from collapsing 
matter. In fact they are solutions of the source-free Einstein or coupled Einstein- 
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Maxwell equations and are characterised for an outside observer completely in 
terms of their mass, their angular momentum and their electro-magnetic charge. 
One might alternatively think of them as black holes that have been formed in the 
infinite past and settled in a stationary state in which all trace of their previous 
history has been lost, e.g. by emitting gravitational radiation. 

Because of this restriction these lectures do not include a detailed account 
of the process of gravitational collapse as studied by astrophysicists, but they 
do allow one to isolate and study in detail the special aspects of the physics of 
black-holes. It may even be that these very special simple solutions are relevant to 
elementary particle physics near the Planck scale, but that leads one to consider 
the problems of quantum gravity, a subject we will not deal with at length in 
these lectures. The interested reader can find more information on many aspects 
of black-hole structure and formation in the literature [28, 29]. 


4.2 The Schwarzschild solution 


The simplest of all stationary black-hole solutions of the source-free Einstein 
equations is that for the static spherically symmetric space-time, asymptotically 
flat at spatial infinity, first described in the literature by Schwarzschild’ [30]. We 
present a derivation of this solution making full use of the symmetries and their 
connection to Killing vectors explained in chapter 2. 

In view of the spherical symmetry we introduce polar coordinates (1,4, y) 
in three-dimensional space and add a time coordinate t, measuring asymptotic 
Minkowski time at r + oo. The Schwarzschild solution is obtained by requiring 
that there exists at least one coordinate system parametrized like this in which 
the following two conditions hold: 

(i) the metric components are t-independent; 

(ii) the line-element is invariant under the three-dimensional rotation group 
SO(3) acting on three-vectors r in the standard (linear) way. 

The first requirement is of course equivalent to invariance of the metric under 
time-shifts t — t + ôt. Thus formulated both conditions take the form of a 
symmetry requirement. 

In chapter (2) we discussed how symmetries of the metric are related to Killing 
vectors. Rotations of three-vectors are generated by the differential operators L 
of orbital angular momentum: 


o o 
Lı = -sine zg — cot 0 cose g L = cos? zg — cot Osiny =, 


(4.1) 
a 


Oy 


'The solution was found independently at the same time in [31]. 


La = 
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Similarly, time-shifts are generated by the operator 


oo 

Ot 
According to the discussion in sect.(2.4) the coefficients of the symmetry opera- 
tors La = L}, are to be components of a set of Killing vectors L4, A = 0,1,2,3: 


Lo (4.2) 


D Las + D, Lay = 0. (4.3) 


These can equivalently be characterised as a set of constants of geodesic flow as 
in eq.(2.45): 


Ja = Lf py. (4.4) 


The relation with the metric is, that the Poisson brackets of these constants of 
geodesic flow with the hamiltonian are to vanish: 


1 
{Ja, H} = 0, H = 5 I” PuPr (4.5) 


The general solution of these conditions for the hamiltonian H is t-independent, 
whilst it can depend on the angular coordinates only through the Casimir invari- 
ant 


L = p 4.6 
Po + sin? 8 ( ) 

It follows that the metric must be of the form? 
ds? = —h? (r)dt? + g'(r)dr? + k?(r)dQ?, (4.7) 


where dQ? = d6? + sin? 0 dy? is the angular distance on the unit sphere, and 
the coefficients h(r), g(r) and k(r)) are functions of the radial coordinate r only. 
Indeed, the hamiltonian derived from this metric is 


1 1 1 1 
H = -— pi — p — ( 2, + 2) 
Dp2 Pt + 2g? Pr + ape P’ + sin? 0 
(4.8) 
L? 
2k?” 


l 4 4 l 3 4 

D2 Pt 2q? Pr 
Note that the set of functions (h,g,k) can not be unique, as one can always 
perform a coordinate transformation r — r(r) to another metric of the same 


general form (4.7) with coefficients (h, g, k). There are several standard options 
to remove this freedom; to begin with, we choose one in which the spherical 


?In this chapter we employ natural units with e = 1. 
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symmetry is manifest and the spatial geometry is as close to that of flat space as 
possible. That is, we take r as the solution of 


k(r) = r’°g’ (r). (4.9) 


Then the line element becomes 


ds? = —W(r)d? + g(r) (dr? + Pa) = (rd? + g°(r)dr?. (4.10) 


Clearly it is rotation invariant, and the three-dimensional space-like part of the 
line element is conformal to that of flat space in polar coordinates, with the 
conformal factor g(r) depending only on the radius, and not on the orientation 
or the time. In the literature the coordinates chosen here to parametrize the 
Schwarzschild metric are called isotropic coordinates. 

Up to the freedom of radial reparametrizations, the form of the line element 
(4.10) is completely determined by the symmetry requirements. To fix the radial 
dependence of the remaining coefficients h(r) and g(r) we finally substitute the 
metric into the Einstein equations. For the Einstein tensor Gy, = Ry, — 1/2 9u R 
we find that it is diagonal with components 


h 2 hg’? 
Gu = 73 (ngr + 2o- +), 
r g 
1 hg? 2 
m = [2g = (hg) |, 4.11 
Ge = = (wa + E + Zay) (11) 
Goe 0 


Goo = sin? ð gh r 

In the absence of sources or a cosmological constant these expressions have to van- 
ish. As boundary conditions at infinity, where space-time must become asymp- 
totically flat, we require that r — oo implies h > 1 and g > 1. Then we have 
three differential equations for two unknown functions, and therefore the equa- 
tions must be degenerate. One can check, that this is the case if the following 
two relations are satisfied: 


hg = g +rg', [r°(hg)'] = 0. (4.12) 
Indeed, in this case the various components of Gy become proportional: 
g 1 
h? r2sin? h “e 


The conditions (4.12) are met for any g(r) of the form 


1 
Grp = Gu = pe Goo = (4.13) 
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b c g+rg 
r) =a+- + — h(r) = ——. 

gira StS. n= È 

It turns out that eqs.(4.14) indeed provide solutions of the Einstein equations 

satisfying the correct boundary conditions. They are most easily obtained from 

the first Einstein equation (4.11): for hg 4 0 it is possible to factor out all 

h-dependence and obtain an equation for g only: 


(4.14) 


2 1 
gg” + “99 = xf = 0. (4.15) 


There is a unique solution with the required normalization for r — oo: 


g(r) = (1 Zy, (4.16) 


where m is an undetermined constant of integration. Eq.(4.12) then immediately 
provides us with the solution for A: 


2r—m 
h(r) = . 
2r +m 
The complete solution for the static, spherically symmetric space-time in isotropic 
coordinates therefore is 


(4.17) 


or — 2 4 
ds? = — E = =) d? + (1 4 z) (dr? + dQ?) (4.18) 
r m r 


Clearly this solution matches Minkowski space at r — oo by construction, and h 
and g change monotonically (decreasing and increasing, respectively) as a func- 
tion of r untill r = m/2. At that value h = 0, and the metric has a zero mode. 
Then the derivation of the solution we have presented breaks down, and its con- 
tinuation to values r < m/2 requires a careful interpretation, discussed below. 
In any case, it turns out that this singularity of the metric is not physical, as 
invariants of the curvature remain finite and well-defined there. Indeed, in the 
following we construct coordinate systems in which the metric is perfectly regular 
at these points and can be continued to regions inside the surface r = m/2. In 
the isotropic coordinate system the apparent singularity of the metric for r = 0 is 
another coordinate artifact. It disappears if one brings spatial infinity to a finite 
distance by an appropriate coordinate transformation. Details are given in the 
next section. 


4.3 Discussion 


Modulo radial reparametrizations r — r(r) the line-element (4.18) is the most 
general static and spherically symmetric solution of the source-free Einstein equa- 
tions matching flat Minkowski space-time at spatial infinity (r — oo). Therefore 
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in the non-relativistic weak-field range (v < 1,r > m) it should reproduce 
Newton’s law for the gravitational field of a spherically symmetric point mass of 
magnitude M, according to which the acceleration of a test particle at distance 
ris 


Ër GM 
TE = — 3 r (4.19) 


As the spherically symmetric solutions of the Einstein equations (4.18) are dis- 
tinguished only by a parameter m, in the Newtonian limit this parameter must 
be related to the mass M, thereby representing the equivalent gravitational mass 
of the object described by the spherically symmetric metric (4.18). 

We first show that our solutions reproduce Newtons law in the limit r >> m, 
v <1 (=c). The equation of motion for a test particle is the geodesic equation 
(2.7). For the spatial coordinates x' this becomes 


Ër’ dx®\* i 
dr? ` (=) Too + (terms ~ O(v/c)) 
(4.20) 
h m ;h 


Note that in the limit v < 1 one has dz? = dt & dr, whilst for r — oo the factor 
h/9? — 1. Upon the identification m = GM we thus indeed reproduce Newton’s 
equation (4.19). 

Taking into account the dimensions of Newton’s constant, as follows from 
(4.19) and is expressed in eq.(1.2), and our natural units in which c = 1, we find 
that m has the dimension of a length, consistent with the observation that m/r 
should be dimensionless. We have already noted that in isotropic coordinates the 
metric component h vanishes at half this distance: r = ro = m/2. Presently this 
only signals that the isotropic system cannot be extended beyond this range, but 
the surface r = ro does have a special physical significance: it defines the location 
of the horizon of a spherically symmetric black hole. This characteristic value 
ro = m/2 of the radial coordinate is known as the Schwarzschild radius. 


For the range m/2 < r < oo the Schwarzschild solution in the form (4.18) is 
well-defined, and it matches the asymptotic Minkowski metric for large values of 
r. We noted that on the Schwarzschild sphere r = m/2 the metric is singular, 
but as a solution of the Einstein equations it can be continued to smaller values 
of r and then it is again well-defined in the domain 0 < r < m/2. However, it 
turns out that this domain of r-values does not in any sense represent the physical 
interior of the horizon, as one might expect naively. Rather, it easy to establish 
that the solution (4.18) is invariant under the transformation 
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m? 


r — J” (4.21) 
which maps the domain (0, m/2) to (m/2,00). Therefore the isotropic coordinate 
system actually represents a double cover of the region between the Schwarzschild 
radius and the asymptotic region r — oo. Note that under this transformation 
h(r) — —h(r), which can be compensated by an accompanying time-reversal 
t — —t. From these observations it follows that the apparent singularity for 
r = 0 disappears as it simply describes the another asymptotic region like that 
for r + oo in the original coordinate system. 

The existence of such a coordinate system presenting a double cover of the 
space outside the black hole might lead one to suspect, that a test particle thrown 
in radially towards the spherical Schwarzschild surface after a long enough time 
would return to an asymptotic Minkowski region. However, this does not happen; 
in the first place, as a simple calculation to be presented later shows, in terms 
of the time measured by an observer at asymptotic infinity it takes an infinite 
period of time for a test particle starting at any finite radial distance r > m/2 
with any finite velocity even to reach the Schwarzschild surface; certainly such 
an observer will never see it start up and set off towards infinity again. Apart 
from this, leaving the Schwarzschild surface would also take an infinite amount of 
asymptotic Minkowski-time. Clearly, these effects are due to the vanishing of h(r) 
on the horizon, which implies that an infinite amount of asymptotic Minkowski 
time t has to pass during any finite period of proper time r. But this is only 
the minor part of the argument; more important is the second reason, to wit 
that time-like geodesics do not flow from values r > m/2 to values r < m/2. 
Instead, as will become clear they flow into a new region of space-time described 
by complex values of r. 

Some calculational details are given later, but the above arguments at least 
indicate that an observer at spatial infinity is essentially disconnected from the do- 
main of space-time beyond the Schwarzschild sphere. Therefore the Schwarzschild 
surface is called a horizon: from spatial infinity one cannot look beyond it. In 
contrast, things are radically different for an observer moving with the test par- 
ticle towards the horizon. For such an observer in free fall, co-moving with the 
test particle on a radial geodesic and starting from an arbitrary point at a fi- 
nite distance from the Schwarzschild surface, only a finite amount of proper time 
passes until the radius r = m/2 is reached, although for large radial distances 
the required time grows linearly with the distance. 

Moreover, having reached the horizon the particle (and the observer) can 
continue falling inwards without noticing anything particular, at least initially. 
The only remarkable effect they could discover once they have passed beyond 
the Schwarzschild surface is the impossibility to turn back to spatial infinity, 
no matter how powerful the engines they have at their disposal to accelerate 
and propel them. Therefore their passage to the inside of the Schwarzschild 
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sphere represents an irreversible event: they are permanently trapped inside the 
Schwarzschild sphere. For this reason the Schwarzschild space-time is called a 
black hole: nothing, not even radiation, can get out once it has fallen through 
the horizon, at least not in the domain of classical relativity. 

To describe this motion of test particles (i.e. the flow of geodesics) through 
the Schwarzschild sphere, one must pass to a different coordinate system, one 
better suited to particles in free fall rather than to the description of the space- 
time from the point of view of an asymptotic Minkowskian observer at rest at 
spatial infinity. We make this passage in two steps. First we show that in other 
coordinate systems there is a different continuation of space-time beyond the 
Schwarzschild surface, one in which the interior geometry inside the horizon is 
distinct from the outside geometry and can not be mapped back isometrically to 
the exterior. After that we show by explicit calculation of the geodesics, that this 
alternative continuation can be taken to decribe the true physical situation for a 
test particle falling through the Schwarzschild sphere. 


4.4 The interior of the Schwarzschild sphere 


To go beyond the region of space-time covered by the isotropic coordinates, we 
perform a radial coordinate transformation which instead of making the spatial 
part of the metric conformal to flat space, as in (4.9), (4.10), makes the angular 
part directly isomorphic to the two-sphere of radius r: 


r = rg(r), k'r) =r. (4.22) 
Then the metric takes the form 


ds? = —h?(r)dt? + g°(r) dr? + Pa’. (4.23) 
The explicit form of the radial transformation (4.22) gives for the coefficients 
- 1 2m 
P = —— =l- L. 4.24 
= FA r (4-24) 


This leads to the standard form of the Schwarzschild solution as presented in the 
original papers [30, 31]: 


2 
ds? = — (1 — m) dt? + dr? + Pd. (4.25) 


1 
eS) 
r 
As to lowest order the transformation (4.22) does not include a rescaling, we have 
for large r and r: 


si| = 
— 
— 


= - + 0 (=). (4.26) 
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Therefore in the asymptotic region the Newtonian equation of motion for a test- 
particle in the field of a point mass M still holds in the same form after replacing 
r by r, and in the non-relativistic long-distance limit the two coordinate systems 
are indistinguishable. 

Note that the Schwarzschild radius r = m/2 is now located at r = 2m. Again 
the metric components are singular there, but now there is a zero mode in the 
(tt)-component and an infinity in the (rr)-component, in such a way that the 
determinant of the metric remains finite on the horizon. On the other hand, 
the determinant g does have a real zero for r = 0. This metric singularity was 
not present in the isotropic coordinate system, as it did not cover this region of 
space-time, but it turns out to represent a real physical singularity: the curvature 
becomes infinite there in a coordinate-independent way. 

A peculiarity of the metric (4.23) is that the components h? and g? are not 
positive definite, but change sign at the Schwarzschild radius. This implies that 
the role of time and radial distance are interchanged for r < 2m: r is the time 
coordinate and t the radial coordinate. At the same time the metric components 
have become explicitly time dependent (they are functions of 7), and therefore 
the metric is no longer static inside the horizon. In return, the metric is now 
invariant under radial shifts t > t+ ôt. Hence it is the radial momentum p; 
which now commutes with the Hamiltonian (4.8), the generator of proper-time 
translations, and is conserved. 

This strange result is the price one has to pay for continuing the Schwarzschild 
geometry to the interior of the Schwarzschild sphere. It is quite clear that this 
could never heve been achieved in the isotropic coordinate system (4.18), at least 
not for real values of r. Together with the appearance of the curvature singularity 
at r = 0 this shows beyond doubt that the region inside the Schwarzschild sphere 
r < 2m has no counter part in the standard isotropic coordinate system. 


To gain a better understanding of the difference between the continuation to 
r < m/2 in isotropic coordinates and r < 2m in Schwarzschild coordinates it is 
useful to perform yet another radial coordinate transformation to a dimensionless 
coordinate p defined by [32] 


tanh É = A(r), (4.27) 
with the result that 
r = rg(r) = 2m cosh? È. (4.28) 
Then the Schwarzschild metric becomes 


ds? = — tanh? E dt? + 4m? cosh“ È (dp? + do?) . (4.29) 
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As cosh p/2 > 1 this expression for the metric, like the isotropic coordinate sys- 
tem, is valid only outside the horizon r > 2m. Moreover, continuing to negative 
values of p we reobtain the double cover of the exterior of the Schwarzschild black 


hole: 


— 2 
tanh (£) = —h(r) = (im), (4.30) 
2 Ar 
Now the interior of the Schwarzschild sphere can also be reparametrized in a 
similar way, taking account of the fact that h? has become negative: 


r = m(1+coso) = 2m cos” = tan? Z = =h? (F). (4.31) 


This parametrization clearly exists only for 0 < r < 2m. In these coordinates 
the metric takes the form 
2 27 42 24T 2 2 

ds“ = tan z% + 4m’ cos 5 (do + dQ ). (4.32) 
The time-like nature of o now is obvious. It is also clear, that this coordinate 
system presents an infinitely repeated periodic covering of the interior of the 
Schwarzschild sphere: the fundamental domain may be chosen to be 0 <0 <r, 
corresponding to 2m > r > 0; but continuation of o beyond this region leads to 
repeated covering of the same set of r-values. 

Comparing these two parametrizations of the exterior and the interior of the 
Schwarzschild sphere, it is now obvious that they are related by analytic contin- 
uation [32]: the interior of the sphere is not described by negative values of p, 
which correspond to r < m/2 in isotropic coordinates, but by imaginary values 
of p: 


o = tip, (4.33) 


which is like a Wick rotation changing the space-like radial coordinate p to the 
time-like coordinate o. In terms of the original isotropic coordinates it can be 
described as the analytic continuation of r to complex values 

r= e®, (4.34) 

2 

Thus in the interior of the Schwarzschild surface the modulus of r is constant, but 
its phase changes, whilst outside its phase is constant and the modulus changes. 
Note that whereas the period of the fundamental domain of ø in describing the 
interior region of the Schwarzschild sphere is 7, the values of the isotropic radial 
coordinate r are periodic in o with period 27. Again it seems that the isotropic 
coordinates cover the Schwarzschild solution twice. We have more to say about 
that in the following. Finally we observe, that as indicated in eq.(4.33) there are 
two ways to perform the analytic continuation from p to ø, 1.e. from the exterior 
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to the interior of the Schwarzschild space-time. Also this fact has non-trivial 
significance for the geometry of the Schwarzschild space-time. 


4.5 Geodesics 


The analysis presented so far makes it clear that the geometry of the Schwarzschild 
space-time is quite intricate. Indeed, at the Schwarzschild sphere r = m/2 (or 
r = 2m and o = p = 0, respectively) several branches of the space-time meet: 
two branches of the exterior connected to an asymptotic Minkowski region, and 
two branches of the interior connecting the horizon with the space-time singu- 
larity. It is therefore of importance to understand the flow of geodesics in the 
neighborhood of the horizon, and see which branches test particles moving in 
the Schwarzschild space-time follow. We begin with time-like geodesics, repre- 
senting the motion of test particles of unit mass. For such a particle one has in 
asymptotic Minkowski space: 


2H = g papy = —l. (4.35) 


This relation was already argued from the definition of proper time in sect.(2.1), 
eq.(2.8) and below. As H is conserved, the relation holds everywhere on the 
particle’s worldline. 

The solution of the geodesic equations for Schwarzschild space-time is much 
simplified by the large number of its symmetries. They provide us with sufficiently 
many additional constants of motion to allow a complete solution of the equations 
of motion of test particles in terms of their energy and angular momentum. The 
first of these constants of motion is the momentum component p;, conserved 
because the metric is t-independent: 


dt dt € 

pi = =h (r) =-« © Zo Rr) (4.36) 
At space-like infinity k? — 1; also, with 7 defined as proper time, dt = dr for a 
particle at rest there. Therefore a particle at rest at r — oo has € = 1, equal to the 
(unit) rest mass of the test particle; moreover, as is intuitively obvious and argued 
more precisely below, if it is not at rest it must have € > 1. Indeed, particles with 
€ < 1 cannot reach spatial infinity and in this sense they are in bound states. 
The factor e gives the usual special relativistic (kinematical) time-dilation for a 
moving particle in asymptotic Minkowski space. 

In contrast, for finite radial distance ro < r < oo the factor h?(r) describes the 
gravitational redshift: the time dilation due only to the presence of a gravitational 
field. It represents a universal effect, not influenced by the state of motion of the 
particle. 

The conservation of angular momentum gives us three conservation laws, 
which are however not independent: 
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J, = —siny pe — coth cosy py, 
Jı = cosy ps — cot sing py, (4.37) 
J3 = Po, 


with J? + J? + J} = L’? =, a constant. Here the separate momentum compo- 
nents are 


dé 
Pe = k?(r) dt’ 


If one orients the coordinate system such that the angular momentum is in the 
x°-direction, then 9 = 7/2 = constant, and therefore pg = 0 and py = J3 = 4. 
Then the angular coordinates as a function of proper time are given by 


. dy 
Po = k’(r)sin? 8 T (4.38) 


d l 
“= (4.39) 


0 = 7/2, dr kr) 


It remains to solve for the radial coordinate? r. For this we use the conservation 
of the hamiltonian, setting 2H = —1: 


dr ? e (2 
27 = |>] =, (—-—- 1). 4.4 
p, (0 a g (i E ) (4.40) 


Note, that consistency of this equation in the region outside the horizon requires 


e L 

2 > Tz +1. (4.41) 
It shows, why in the region r — oo the time dilation factor has to be greater 
or equal to unity: € > 1. At the same time, for non-zero values of angular 
momentum f the radial distance r can in general not become arbitrary low: like 
in the ordinary Kepler problem there is a centrifugal barrier to be overcome. For 
details we refer to the literature [11, 28, 29]. 

Combining eq.(4.40) with that for the angle y, by eliminating r we directly 

obtain the orbital equation 


OE fe Eg (4.42) 


It is clear that an orbit can have a point of closest approach, or furthest distance, 
only if dr/dy = dr/dr = 0. The only physical solutions are points where relation 
(4.41) becomes an equality. When the equality holds identically, for all points of 
the orbit, the orbit is circular. One can show that the smallest radius allowed for a 


3Note that each of these equations is still true for arbitrary radial parametrization. 
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stable circular orbit corresponds in Schwarzschild coordinates to r = 2ml? = 6m. 


To the set of time-like trajectories we add the light-like geodesics, representing 
the orbits of massless particles. In that case the proper-time interval vanishes 
identically on any geodesic, that is to say H = 0. We parametrize the trajectories 
with an affine parameter À and obtain 


dt 
= -k — = — 4.43 
allowing one to replace the affine parameter by t through 
h2 
dA = — dt. (4.44) 
Y 


One can again choose the plane of motion to be 0 = 7/2, and the conservation 
of angular momentum holds in the form 


2 
yi de _ 2de 
L Z = pR 4.45 
hd | dy (4-45) 
where w is a constant. The radial equation now becomes 
dr\* do\* dy\? 
2 2 7 24K 2.52 a 
wae (F) +e (F) + esto (E) (4.46) 
or using the choice 0 = 7/2 above: 
dr\* wh? k? 
*(—) = k (1 -——_}. 4.4 
É (z) ( yr aa) 


Like in the massive case, if w Æ 0 one can eliminate the affine parameter À also 
in favour of y instead of t. 


An interesting case, as concerns the exploration of the horizon and the connection 
between the interior and the exterior of the Schwarzschild sphere, is that of radial 
motion: l = 0 or w = 0, respectively. Let us first consider the time-like geodesic 
for a particle starting from rest near r = oo, which implies € = 1 and dr/dt < 0. 


Then 


dt 1 dr 1 

— = — — = —-— V1 —h?. 4.48 

dT h2’ dT gh ( ) 
It is easy to solve these equations in the Schwarzschild coordinates; indeed, with 
(1 — h?) = 2m/r one finds a solution valid in the whole range 0 < F < co, both 
outside and inside the horizon: 
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2 
3 
A) _ [—(-7)] (4.49) 
where To is the proper time at which the particle reaches the singularity at r = 0. 
Taking this time as fixed, r denotes the proper time, prior to 7, at which the 
particle is at radial distance r. Observe, that this time is finite for any finite 
radial distance r; only when starting from infinity it takes an infinite proper time 
to reach the singularity. Nothing special happens at the horizon r = 2m: it is 
crossed at proper time 


Am 
3? 
a proper period 4m/3 before the singularity is reached (of course, in terms of the 
time-coordinate 7 inside the horizon, it takes a finite period 2m). 

Solving for t as a function of proper time, we have to distinguish two cases: 


(4.50) 


TH = To — 


r > 2m andr < 2m. Namely, introducing the new variable 


y = P (To — ae (4.51) 


we obtain the differential equation 
—- — — = —— 4.52 
4m dr y2- 1’ ( ) 


the solution of which depends on whether the r.h.s. is positive or negative. We 
find for the two cases the solutions 


wje 


i 3 
to —t Tt) —T | 3 1 arccoth = (To — r)| , r>2m; 


w]e 


3 
arctanh + (To — r)| , r <2m. 
4m 
(4.53) 
In both cases the left-hand side becomes infinite at the horizon, i.e. at proper 
time Tg. This confirms that an external observer at infinity will never actually 
see the particle reach the horizon. Similarly, the distance measured from the 
horizon to the singularity at r = 0 is infinite in terms of the coordinate t. 
Now consider what happens in isotropic coordinates. In terms of a new vari- 


able 
2r 
=,/— 4.54 
z =Z, (4.54) 


where we choose the positive branch of the square root, integration of the equation 
for r(7) is elementary and gives 


67 


2 1 1\3 
- (m—-T) = 3 (« + -) . (4.55) 
Clearly this is symmetric under interchange x — 1/zx, corresponding to inversion 
w.r.t. the horizon r = m/2. Indeed, both at r — oo and r = 0 we see that 
T — —oo. This confirms our earlier remark that a particle released from space- 
like infinity does not return there. Instead, it will cross the horizon z = 1 at 
T = Ty, as in the previous calculation. Then continuing r to complex values 


(4.34), we find 


r = e? , = — cos’ —. (4.56) 


The singularity o = a is then reached for T = To, as expected. In terms of the 
radial coordinate p, related by the pseudo-Wick rotation (4.33) the equation for 
the radial distance in terms of proper time outside the horizon becomes 

T —T 1 


m, 3 cosh? E, (4.57) 


The analytic continuation of coordinates in the complex plane on crossing the 


horizon is manifest. 
Next we consider incoming light-like radial geodesics. The equation is 


— = —-. 4.58 
Let to refer to the time at which the massless particle (e.g., a photon) passes a 
fixed point on the radius, e.g. ro in Schwarzschild coordinates. Then the solution 
of the geodesic equation for t is 

t—to oF 


= — In 


2m 2m 


(4.59) 


With the absolute value as argument of the logarithm this equation is valid both 
outside and inside the horizon. The time for the photon to reach the horizon 
from any point outside diverges, as does the distance it has to move from horizon 
to the curvature singularity, as measured by the coordinate t. 


4.6 Extended Schwarzschild geometry 


As we have already noticed that for a massive particle falling across the horizon 
nothing very special happens, it is unsatisfactory that we do not have a non- 
singular description of the same process for massless particles. To improve the 
description of light-like geodesics at the end of the last section, it is not suffi- 
cient to perform a radial coordinate transformation: one has to eliminate the 
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t-coordinate and replace it by one which does not become infinite on the horizon. 
Moreover, to describe light-like geodesics it would be convenient to have a metric 
where the light-cone structure is manifest by having the radial part of the met- 
ric conformal to radial Minkowski space: a coordinate system in which (r,t) are 
replaced by coordinates (u,v) such that the radial part of the metric becomes 


—h?(r)dt? + g(r} dr? + f(u,v) (-dw? + dv”) . (4.60) 


Such a coordinate system was first proposed by Kruskal and Szekeres [33]; defining 
the new coordinates such as to be dimensionless, the line element is written in 
the form 


ds? 


4m? 


= F(u,v) ( dv’ 4 du?) h g(u,v) dQ’, (4.61) 


with the standard choice for (u,v) defined in terms of the Schwarzschild coordi- 
nates (r,t) by: 


t 
tanh —, if|v] < |u|; 
2 2 r = u 4m 
u — v = | — — 1) em, — = (4.62) 
v t 
coth —, if|v] > lul. 
Am 


In these coordinates the horizon r = 2m, t = too is located at u = +v, whilst 
the curvature singularity r = 0 corresponds to the two branches of the hyperbola 
u? — v? = —1, one in the past and one in the future of an observer at space-like 
infinity. More generally, from the first result one concludes that the surfaces r = 
constant are mapped to hyperbola’s u?—v? = constant, where the last constant is 
positive outside and negative inside the horizon. These hyperbola’s are worldlines 
of particles in circular orbit, i.e. being subject to a constant acceleration. 

Similarly hypersurfaces t = constant correspond to u = constant x v, repre- 
sented by straight lines in a (u,v)-diagram. Note also that the singularity r = 0 
consists of a set of two space-like surfaces in (u,v)-coordinates. 


More explicitly, we can write (u,v) in terms of (r,t) as 


r s t 
u = +),/— — lem cosh —, 
2m Am 
r r t 
v = 4+,/——lem sinh —, 
2m Am 


outside the horizon (F > 2m); and 
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(4.64) 


To r t 
v = +,/1 ——e™ cosh —, 
2m Am 


within the horizon (r < 2m). We note once again, that these coordinates produce 
two solutions for each region of the space-time geometry. In the above conventions 
the dimensionless coefficient functions are given implicitly by 


flu, v) = > © am glu, v) = Am2 (4.65) 


We can also compare with the coordinate systems (4.29) and (4.32), describing 
explicitly the double cover of the Schwarzschild space-time: 


t t 
u = er hE sinh É cosh —, v = 62 oF sinh? sinh —, (4.66) 
4m 2 4m 


for the exterior region |v| < |u|; and 
t 


. oO , 
sin — sinh —, 
2 Am 


on 


t 
sin 5 cosh In’ (4.67) 


2 


1 z 2 
z3 COs 2 


1 fom 
2 cos 2 


for the interior region |v| > |u|. From these expressions we infer, that the exterior 
regions p > 0 and p < 0, corresponding to r > m/2 and r < m/2 in isotropic 
coordinates, and the two interior regions o > 0 and g < 0 as well, coincide with 
the regions (u,v) > 0 and (u,v) < 0 in Kruskal-Szekeres coordinates, respectively. 
Hence the double cover of Schwarzschild space-time given by the two solutions for 
the Kruskal-Szekeres coordinates above correspond exactly to the double cover 
we found in isotropic and (p, 7 )-coordinates. 

The Kruskal-Szekeres coordinates were introduced, and are especially useful, 
because they represent the light-like geodesics in a very simple way. In particular, 
the radial light-like geodesic are given by 


2 2 du 
duu = d? => T ~ +1. (4.68) 
These correspond to straight lines in the (u, v)-diagram parallel to the diagonals. 
The special light rays u = +v are the asymptotes of the singularity and describe 
the horizon itself: they can be interpreted as photons permanently hovering on 
the horizon. Any other light-like radial geodesic eventually hits the singularity, 
either in the past or in the future, at finite (u,v), although this is in the region 


before t = —oo or after t = œ in terms of the time of an observer at spatial 
infinity. 
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As the light cones have a particularly simple representation in the (u,v)- 
diagram, for all regions inside and outside the horizon, it follows that time-like 
and space-like geodesics are distinguished by having their derivatives |du/dv| > 1 
for the time-like case, and |du/dv| < 1 for the space-like ones. One can check 
by explicit computation, that a particle starting from rest at any distance r and 
falling into the black hole has a time-like world line everywhere, also inside the 
horizon: du/dv can change on the world line of a particle, but never from time- 
like to space-like or vice versa. Therefore all time-like geodesics must cross the 
horizon (light-like) and the singularity (space-like) sooner or later. If they don’t 
reach or start from spatial infinity, they must come from the past singularity 
across the past horizon, and then cross the future horizon to travel towards the 
final singularity, all within a finite amount of proper time. 

The Kruskal-Szekeres coordinate system also teaches us more about the topol- 
ogy of the Schwarzschild geometry: how the various regions of the space-time are 
connected. First of all it shows that there is a past and a future singularity inside 
a past and a future horizon (|u| < |v|), depending on v being positive or negative. 
In this respect it is helpful to observe, that according to eq.(4.64) for 7 < 2m the 
sign of v is fixed for all t. In the literature these regions are sometimes refered to 
as the black hole, the interior of the horizon containting the future singularity, and 
the white hole containing the past singularity. In a standard convention, these 
regions of space-time are classified as regions I (future Schwarzschild sphere) and 
IV (past Schwarzschild sphere), respectively. 

There are also two exterior regions connected to asymptotic Minkowski space. 
They are distinguished similarly by u > 0 and u < 0, as for r > 2m the sign of u is 
fixed for all t: one cannot reach negative u starting from positive u or vice versa 
by any time-like or light-like geodesic. In the same standard convention these 
regions are refered to as region I and region III, respectively. Being connected by 
space-like geodesics only, the two sheets of the exterior of the Schwarzschild space- 
time are not in causal contact and cannot physically communicate. However, the 
surfaces t = constant do extend from u = œ to u = —oo; indeed, the regions 
are space-like connected at the locus u = v = 0, where the two branches of the 
horizon intersect. To understand the topology of Schwarzschild space-time, we 
must analyze this connection in some more detail; this can be done using the 
technique of embedding surfaces, as explained e.g. in ref.[11]. 

Consider a surface t = constant and fix a plane 0 = 7/2. We wish to un- 
derstand how the regions u > 0 and u < 0 are connected for this plane. The 
Schwarzschild line element restricted to this plane is 


ds? = om, + rdg’. (4.69) 


Now we embed this two-dimensional surface in a three-dimensional flat, euclidean 
space: 
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ds? = dz? + dy? + dz’, (4.70) 
by taking z =rcosy, y = Tsin y, and defining the surface z = F (r) such that 
ds? = (14 F?(r)) dr? + Pde. (4.71) 
Comparison of eqs. (4.71) and (4.69) then implies 


2 o) -> ı, (4.72) 


16m? 4m 2m 


This defines a surface of revolution, obtained by taking the parabola z? = 8mg — 
16m? in the (x, z)-plane and rotating it around the z-axis. It is a single surface, 
with looks like a throat, of which the upper half z > 0 corresponds to 


5,071 29. or uJ, (4.73) 


whilst the lower half z < 0 describes the branch 


5-7! <0, or u<0, (4.74) 


The two regions are connected along the circle r = 2m at z = 0. Of course, the 
spherical symmetry implies that the same topology would be observed for any 
choice of angle 0; any plane with u > 0 is connected to one with u < 0 in this way. 
Although no time-like travel is possible classically between the regions on each 
side of the throat, the possibility of particles tunneling via quantum processes 
between the two regions poses a curious and interesting question to theories of 
quantum gravity. 


4.7 Charged black holes 


The spherically symmetric solution of the source-free Einstein equations can be 
extended to a solution of the coupled Einstein-Maxwell equations including a 
static spherically symmetric electric field (in the rest frame of the black hole)[34]. 
As the only solution of this type in flat Minkowski space is the Coulomb field, we 
expect a solution of the coupled equations to approach the Coulomb potential at 
space-like infinity: 


q 
A = A ,dz" = dt 4. 
H T 4rrp(r) 9 ( 75) 


with p(r) > 1 for r + co. This is a solution of the Maxwell equation 


D,F™ = 0, (4.76) 
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in a spherically symmetric static space-time geometry. The Maxwell field strength 
appearing in this equations is 
q prp 


1 y p 
F = z Pw da" A de = 2p? dt A dr. (4.77) 


In isotropic coordinates (4.10) the energy-momentum tensor of this electro-magnetic 


field reads 


2 2 


1 ¢ (prp) ~ fh . 
— A 2) _ 2 2 orn? 
TalF] = (FaF, — F Juv”) = 32 app tae pr „r“ sinf0 |. 
(4.78) 
This has the same general structure as the Einstein tensor (4.13): 
1 hg? 2 h? 
Gu = oh (ws + 7 += (oa!) diag (5 —1,r?,r? sin? s) , (4.79) 


which was derived under the conditions (4.12), with the general solution (4.14). 
Now requiring the coupled Einstein-Maxwell equations 


Gu = -81G T [F], (4.80) 


this reduces to a single condition on the radial functions (h, g, p), with the solution 


m V 
pr) = gir) 1 2 
r Ar 
(4.81) 
gt+rg Ar? — p? 
h(r) = = ——__. 
g 4r? + Amr + v7? 
Here the constant v? is defined by 
2G 
v= m—e’, e = ile (4.82) 


Ar 


It may be checked that this solution satisfies the Maxwell equations in this grav- 
itational fields as well. As g(r) can also be written in the form 


m\? e? 
={1l4+—] - — 4. 
gir) = (+2) - 5, (4.83) 
we find that the Schwarzschild solution is reobtained in the limit e? — 0. For 
v? > 0 the zero of the metric coefficient k? is now shifted to 
v 
= —. 4.84 
r= (4.84) 
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As in the Schwarzschild case this is a point where two asymptotically flat regions 
of space-time meet: the metric is again invariant under the transformation 


py? 


r> g t > -t, (4.85) 
under which h(r)dt + h(r)dt. Moreover, also like in the Schwarzschild geometry, 
the sphere r = v/2 actually defines a horizon and time-like geodesics are contin- 
ued into the interior of this horizon, which is not covered by the (real) isotropic 
coordinate system. 

A Schwarzschild-like coordinate system which allows continuation into the 
interior of the horizon exists. It was found by Reissner and Nordstrom [34] and 
is related to the isotropic coordinates by a radial coordinate transformation 

„2 

rF=rgr)=rtm+t T (4.86) 
r 

Then the metric of the charged spherically symmetric black hole takes the form 


d? = —h?dt? + gdr? + r°dQ’, (4.87) 


where the coefficients (h, g) are again inversely related: 


2 2 2 
hE) = — =1-45-(-2) - Z, (4.88) 
r r 
The horizon is now seen to be located at r = m + v; in fact, there is a second 
horizon at F = m — v, provided m > e. As in general h? and g? are not posi- 
tive definite, at each of these horizons the role of radial and time coordinate is 
interchanged. Only for v = 0, i.e. for m = e, the situation is different, as the two 
horizons coincide and no interchange of time- and space-like coordinates occurs. 
This special case of Reissner-Nordstrom solution is refered to in the literature as 

an extremal charged black hole. 

The result of the second horizon for non-extremal charged black holes is, 
that the real curvature singularity that is found to exist at r = 0 is a time- 
like singularity, for m > e. This is expected from the Coulomb solution in flat 
Minkowski space, which contains a time-like singularity of the electric field at 
r = 0. Therefore such a black hole looks more like a point particle than the 
Schwarzschild type of solution of the Einstein equations. 

The interior region between the horizons can be obtained from the isotropic 
coordinates by complex analytic continuation of r from the outer horizon by 

Vis 

r= je. (4.89) 

The domain of the argument o is chosen as [0,7]. Indeed, for these values of r 
the radial Reissner-Nordstrom coordinate is 
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r = m + -vcoso, (4.90) 


and takes all values between the two horizons r = m+ v precisely once. In this 
parametrization the isotropic Reissner-Nordstrom metric becomes 


v? sin? o 


ds? = dt? + (m + v cos a) (—do? + dO?) . (4.91) 


(m + v cos a) 


Again we observe the role of o as the time-like variable, with t space-like between 
the horizons. For o = m the isotropic coordinate r becomes negative: 


V 
C=T>r= -7 (4.92) 


Thus we discover that the Reissner-Nordstrom space-time can be extended con- 
sistently to negative values of the radial coordinate. In particular the curvature 
singularity at r = 0 is seen to be mapped to 


more 
2 


This shows first of all, that there are two (time-like) branches of the singularity, 
and secondly that r can become negative for non-extremal black holes with m > e. 


r= — 


(4.93) 


The two branches are found also upon analytic continuation of ø in (4.91) to 
imaginary values: 


0 — +ip. (4.94) 


This leads to a parametrization of the Reissner-Nordstrom geometry for the ex- 
terior regions of the form 


v? sinh? p 


ds? = — 


dt? + (m v cosh p} (dp? + dO?) , (4.95) 


(m v cosh p)” 


where the plus sign holds for the exterior region connected to asymptotic flat 
Minkowski space, and the minus sign to the region connected to the curvature 
singularity. Because of the double-valuedness of the analytic continuation (4.94) 
there are two copies of each of these exterior regions involved. 


4.8 Spinning black holes 


The Schwarzschild and Reissner-Nordstrom black holes are static and spheri- 
cally symmetric w.r.t. an asymptotic Minkowski frame at spatial infinity. There 
also exist extensions of these solutions of Einstein and coupled Einstein-Maxwell 
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equations which are neither static nor spherically symmetric, but represent ro- 
tating black holes with or without charge [35, 36]. These solutions have an axial 
symmetry around the axis of rotation, and are stationary: in the exterior region 
the metric is time-independent, and the rate of rotation and all other observable 
properties of the black hole are constant in time. 

The standard choice of coordinates for these metrics, with the axis of rotation 
being taken as the z-axis and reducing to those of Schwarzschild, or Reissner and 
Nordstrom for the spherically symmetric case of no rotation, is that of Boyer 
and Lindquist [37]. Denoting these coordinates by (¢,7,9,y~) to emphasize their 
relation with the spherical coordinate systems for the non-rotating metrics, the 
line element for reads 


ds? = A (at — asin? ody) + 


r2 + a? cos? 6 


r? + a? cos? 6 a x2an2 
— (ar + A*dé ) 

F? + a’)? sin? 0 a 2 

+ U ta) sin’ 6 (up _ —"— it) , 
r2 + a? cos? 0 r2 +a? 
(4.96) 

Here a is a constant parametrizing the deviation of this line element from the 
usual diagonal form of the Schwarzschild-Reissner-Nordstrom metric; its physical 
interpretation is that of the total angular momentum per unit of mass: 


J = ma. (4.97) 


As one can always choose coordinates such that the angular momentum points 
along the positive z-axis, the angular momentum per unit of mass may be taken 
positive: a > 0. For convenience of discussion this will be assumed in the follow- 
ing. 

The quantity A? in the line element (4.96) is short hand for 


A? = (r m)’ La Heem. (4.98) 
Note that, like the metric coefficients h? and g? we have encountered earlier, 
this quantity is not positive definite everywhere. Also, the zero’s of this function 
determine the location of the horizons, of which there are two (an outer and an 
inner one) for m? > a? + e?: 


Thy =moa Vm? — a? — e. (4.99) 


The gravitational field represented by the metric (4.96) is accompanied by an 
electro-magnetic field, the vector potential of which can be chosen as the one- 
form 


qr . 2 
A= Tap (at — asin ode) . (4.100) 
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For convenience of notation we have here introduced the function p defined as 


P = r + acos? ð. (4.101) 
These fields are solutions of the coupled Einstein-Maxwell equations (4.76) and 
(4.80), with the Maxwell field strength given by the two form 


F = Ta (P — a’ cos? 0) dr ^ (at — asin? ode) 


(4.102) 


qa Tcos0sinð 


27 pt 


d9 ^ (adt — (7? + @è)dẹ), 


and with the identification of the quantity e? as in eq.(4.82): 


To determine the properties of the space-time described by the Kerr- Newman line 
element we need to analyse the geodesic flow. Because the spherical symmetry 
of the Schwarzschild and Reissner-Nordstrom solutions is absent, and only axial 
symmetry remains, there is one less independent constant of motion from angular 
momentum: only the component along the axis of rotation (the z-axis, say) is 
preserved. Fortunately this is compensated by the appearance of a new constant 
of motion, quadratic in momenta, which replaces the Casimir invariant of total 
angular momentum [39, 40]. As a result one can still completely solve the geodesic 
equations in terms of first integrals of motion. We now give some details. 

First we write down the geodesic hamiltonian, effectively representing the 
inverse metric g””: 


1 — . 1 2 1 7 2 
H = 37? are + p + (asind + = pe) — Kz ((7? +a?) pi + apo) | 
(4.103) 


Clearly, to be well-defined we must assume 9 # (0,7), as the z-axis represents 
a coordinate singularity in the Boyer-Lindquist system where the angle y is not 
well-defined. The above hamiltonian is a constant of geodesic flow; recall, that for 
time-like geodesics it takes the value 2H = —1, whilst for light-like ones H = 0. 

Next there are two Killing vectors, respresenting invariance of the geometry 
under time-shifts and rotations around the z-axis: 


dt dy 
E€ = -p = TITIA = Ite Dy? 
(4.104) 
dt dy 


C= Po = Jet Ty + Gee Ty 
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Of course, for time-like geodesics the affine parameter à maybe taken to be proper 
time 7. These equations can be viewed as establishing a linear relation between 
between the constants of motion (€, £) and the velocities (dt/d\, dy/dX): 


ORES (4.105) 


with the dot denoting a derivative w.r.t. the affine parameter A. The determinant 
of this linear transformation is 


det G = —gugeo + J, = A? sin? 0, (4.106) 


and changes sign whenever A? does. Therefore the zero’s of A? define the locus 
of points where one of the eigenvalues of G changes sign, from a time-like to 
a spacelike variable or vice-versa; at the same time, the coefficient of p? in 
the hamiltonian (4.103) changes sign, confirming that these values represent the 
horizons of the Kerr-Newman space-time. It follows from this analysis, that 
outside the outer horizon, in the region connected to space-like infinity, det G>0. 

Now if from eqs.(4.104) the angular proper velocity dy/dA is eliminated, one 
finds the relation 


det & dt 
e- MDZ (4.107) 
Jop dà 
where the quantity Q with the dimensions of an angular velocity is defined by 
Q = -H an tt P-A > 
Jop (T? + a?)? — A?a? sin“ 0 
(4.108) 
a 1 


Fpa? 1+4 A2p? (P2 + a2)-! (P + a? ANE 
Observe, that it is proportional to a with a constant of proportionality which 
for large but finite r is positive. Hence 2 represents an angular velocity with 
the same orientation as the angular momentum J. With a canonical choice of 
coordinate system this quantity is positive, and therefore Q > 0 as well. 

It also follows from eq.(4.107) that for time-like geodesics outside the horizon, 
with A = 7 and therefore dt/dà > 0, an inequality holds of the form 


e- >O. (4.109) 


For large r near infinity one finds in fact the stronger inequality € — Q£ > 1. 

Returning to the solution of the geodesic equations, we take into account the 
existence of a constant of motion K such that 2K = K“ p,p, with K*” the 
contravariant components of a Killing tensor: 


Du Koy = 0. (4.110) 
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The explicit expression for K is 


> 1 2.2 2 2 2 r’ sin’ 4 Pe 
K = zz (Ae cos“ 6 p + Fp) + 252 (an+ P55) 
(4.111) 
a? cos? 0 / _ ? 
+ Soar (+e ype + apo) 


That K defines a constant of geodesic flow is most easily established by checking 
that its Poisson bracket with the hamiltonian H (4.103) vanishes: 
{k,H} = 0. (4.112) 


Using these constants of motion, it is straightforward to establish the solutions 
for time-like geodesics in the form 


dt r+ a? (+a F — A?a? sin? 6 
ae _ / LTr Y T’, naam 
P dr a ( A2 tE A2 , 
ple _ pXe a? sin? 8 al +a’ 

dr A? sin? 6 A? í 


(4.113) 
—-A? (2K +7) + (la -e +a°)) 


D; 
A 
[N 
`a 
|S 
OO 
bo 
II 


2 2 

p“ (2) 2K — a? cos? 6 — a’ sin? 0 ( < a -) . 

Note that for a particle starting from r = œœ one has ¢ > 1. These solutions are 
well-defined everywhere except on the horizons and at points where p? = 0. At 
the latter locus the curvature invariant Ravos R'”’? becomes infinite and there is 
a real physical singularity. At the horizons both ¢ and y change infinitely fast, 
signalling infinite growth of these coordinates w.r.t. an observer in asymptotic 
flat Minkowski space. 

As with the infinite amount of coordinate time ¢ to reach the horizon in the 
Schwarzschild and Reissner-Nordstrom solutions, this merely signifies a coordi- 
nate singularity. However, here it is not only t but also y which becomes infinite. 
To gain a better understanding of the reason for this rotation of geodesics, first 
notice that according to the second equation (4.113) even for £ = 0 the proper 
angular velocity dy/dr does not vanish in general, because € 4 0. Hence there are 
no purely radial geodesics: any incoming test particle picks up a rotation. Now 
computing the angular velocity as measured by an observer at rest at spatial 
infinity from the first two equations (4.113) one finds: 
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dy _ Eft T lgu 


= . 4.114 
dt Edo + liio ( ) 


It follows immediately, that for a geodesic without orbital angular momentum: 
L = 0, the observed angular velocity equals the special reference value (4.108): 


Z| =Q. (4.115) 


(4.116) 


Comparison with the expression for Q, eq.(4.108), shows that outside the horizon 
the inequality 


Q < Qp, (4.117) 


holds, with equality only on the horizon. As these quantities are finite, the 
time ¢ can increase near the horizon at an infinite rate only if y increases at a 
proportional rate. For this reason both ¢ and y have to grow without bound as 
r approaches r4p. 

At space-like infinity the Killing vectors (4.104) define time translations and 
plane rotations. However, their nature changes if one approaches the horizon. 
Consider the time translation generated by p; = —e. The corresponding Killing 
vector has contravariant and covariant components, given by 


E =O S Eip = — Oty (4.118) 


Here 6" is the usual Kronecker delta symbol. Similarly the components of the 
Killing vector for plane rotations, generated by pọ, are 


b= 85 bon = Ion (4.119) 
It follows, that these Killing vectors have norms given by 
7 —A? 4+ a? sin? 0 


(&) = ga = z ; (4.120) 


and 


(r? + a7)? — A?2a? sin? 6 


(Eo) = gop = sin” 8 7 , (4.121) 


whilst their inner product is 
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F +a- A? 


532 


p 
From these equations several observations follow. First, as the norms and in- 


(4.122) 


Ert Eo = —Gty = asin? 6 


ner products of Killing vectors are coordinate independent quantities, we can 
write the expressions for the reference angular velocity Q, eq.(4.108), and for the 
observed angular velocity, eq.(4.114), in a coordinate independent form: 


br: Eo do &-&e-~ GE 
Q = , — = 4, 4.123 
E dt  &E— Eril l ) 
The last equality is equivalent to 
de @ @ 
(e —OQl) — = Ne —- (=D fe -0E . (4.124) 
dt & VE 


Second, with a > 0 and excluding the z-axis 9 = (0,7) —where the Boyer- 
Lindquist coordinates are not well-defined and Q vanishes— Q and the invariant 
inner product é -ép are positive outside the horizon. Also the positivity of det G 
in the exterior region, eq.(4.106), can be expressed in the invariant form 


Q? > Si = ‘i <1 (4.125) 
Z @ “YE > t+) 

with the equality holding only on the horizon. Then for £ > 0 we find from 
eq.(4.124) that dy/dt > Q, whilst for £ < 0 the angular velocity is less than this 
reference value: dy/dt < Q. As long as é is time-like this is all that can be 
established; in fact, with l < 0 and €? < 0, the absolute value of the angular 

velocity dp/dt can become arbitrarily large. 
However, it is easily established that the norm of € changes sign on the surface 


r(0) = 7o(8) = m + Vm? — e? — a? cos? 0, (4.126) 


which, with the exception of the points where the surface cuts the z-axis: 9 = 
(0,7), lies outside the horizon. The region between this surface and the horizon 
is called the ergosphere, and é is space-like there: €? > 0. It then follows directly 
from eqs.(4.124) and (4.125) that in this region the angular velocity is always 
positive: 
0 org 

T: > 0, or & > 0. (4.127) 
Thus we have established, that inside the ergosphere any test particle is dragged 
along with the rotation of the black hole and no free-falling observer can be non- 
rotating there. Related to this, both & and €, are space-like Killing vectors there. 
Nevertheless, outside the horizon there still exists time-like combinations of these 
vectors even inside the ergosphere: 
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X = & — Ey, (4.128) 


which has non-positive norm: x° < 0 as a result of the inequality (4.125), with 
equality holding only on the horizon r = ry. 


4.9 The Kerr singularity 


We continue the analysis of black-hole geometry with a brief discussion of the 
nature of the singularity of spinning black holes. To begin, we have already 
noticed that, like the spherical charged black hole, the spinning Kerr-Newman 
black holes have two horizons, an innner and an outer one. As the singularity 
p = 0 is located inside the inner horizon, we expect the singularity to have a 
time-like character. 

In order to elucidate the topological structure of the space-time near the 
singularity, we introduce new coordinates in two steps. First we define the Kerr 


coordinates (V,) as the solution of the differential equations 
(+a) 
dV = dt + AD 


With these definitions the Kerr-Newman line element becomes 


_ ~ a 
dr, dp = dp + Ka T- (4.129) 


. 2 
ds? = (av — asin? odg) + sin 0 


-7 


(P + a’\dp — ad) 
(4.130) 


+2 (dV — asin? 9 dẹ) dr + P d8. 


Now we turn from quasi-spherical coordinates to quasi-Cartesian ones, known as 
Kerr-Schild coordinates, by defining 


x = (Fcosp~—asing) sinð, 
y = (Fsing+acos¢) sind (4.131) 
z = rcos8, 


and we reparametrize the time coordinate t > ¢ using 


di = dV — dr. (4.132) 


We then find that at fixed (0, $) the transformation 
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xr > wv = 2 — 2r cosġ, 
y > y=y - sing, (4.133) 
z >œ 2g = Z, 

is equivalent to r > r’ = —r. This shows, that like the Reissner-Nordstrom case, 


the negative values of r are to be taken as part of the physical domain of values. 
Next one introduces the coefficients k’ (i = 1,2,3) by the relation 


kdr = Oe i, (4.134) 


With these definitions the metric can be written in the hybrid form 


ds? = dz? + u(r) (k- dr)’, (4.136) 


where the first term on the right-hand side is to be interpreted as in flat space- 
time, with x“ = (x,y,z,t). Therefore u and k parametrize the deviation from 
flat space-time. 

Having introduced this new parametrization, the physical singularity p? = 0, 
equivalent to r = 0 and @ = 7/2, is now seen to be mapped to the ring 


r? +y =a, z = 0, (4.137) 


with the solution 


xr = —asing, y = acos¢. (4.138) 


A consequence of this structure is, that in contrast to the Schwarzschild geometry, 
geodesic crossing the horizon do not necessarily encounter the singularity, but can 
pass it by. 


4.10 Black-holes and thermodynamics 


When a test particle falls into the ergosphere of a black hole, the constant of 
motion €;-p = € becomes a momentum component, rather than an energy. 
This provides a means of extracting energy from a spinning black hole, as first 
proposed by Penrose [41]: if a particle enters the ergosphere with energy €, this 
becomes a momentum; during its stay in the ergosphere the momentum can be 
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changed, for example by a decay of the particle into two or more fragments, in 
such a way that the final momentum of one of the fragments is larger than the 
original momentum. If this fragment would again escape, which is allowed as 
the ergosphere is located outside the horizon, it leaves with a different energy 
e’ corresponding to its changed momentum. This energy is now larger than the 
energy € with which the original particle entered. By this process it is possible 
to extract energy from a black hole, even in classical physics. 

It would seem that the Penrose process violates the conservation of energy, 
if there would be no limit to the total amount of energy that can be extracted 
this way. However, it can be established that such a violation is not possible and 
that the process is self-limiting: emission of energy is always accompanied by 
a decrease in the angular momentum (and/or electric charge) of the black hole 
in such a way, that after extraction of a finite amount of energy the black hole 
becomes spherical. After that there is no longer an ergosphere and no energy 
extraction by the Penrose process is possible. 

We will now show how this comes about for the simplest case of an electrically 
neutral black hole of the Kerr type. In this case the location of the horizon is at 


1 
ry = Mm + Vm? — a = m + — vym- J. (4.139) 
m 


Then the angular momentum of test particles moving tangentially to the horizon 
is 


Q=- < = = l (4.140) 
H Rpa 2m m? + /m4t — J? l 
This quantity can be interpreted in terms of the area of the spherical surface 
formed by horizon of the black hole: 


_ Ara 


Ay = Tf dode gages = 4r (Fi, +07) = T (4.141) 
r=rH H 


Eq.(4.140) thus provides an expression for the area of the horizon in terms of the 
mass and the angular momentum of the black hole: 


au = m + Vmi- P. (4.142) 
T 


Now we observe that 
(i) the area Ap increases with increasing mass; 
(ii) at fixed area Ap the mass reaches a minimum for J = 0; 
(iii) by absorbing infalling test particles (including Penrose-type processes) the 
black-hole area Ap never decreases. 

The first two statements are obvious from eq.(4.142). The last one follows if 
we consider the effect of changes in the mass and angular momentum of the black 
hole. We assume that the absorption of a test particle with energy ¢ and angular 
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momentum £ changes a Kerr black hole of total mass m and angular momentum J 
to another Kerr black hole of mass m’ = m+e and angular momentum J’ = J +4. 
This assumption is justified by the uniqueness of the Kerr solution for stationary 
rotating black holes. 

Now under arbitrary changes (dm, óJ) the corresponding change in the area 


SA 2mdm 2 JSJ 
Er (im tm E 


J 
Qyvy mi — J? l 


For an infalling test particle only subject to the gravitational force of the black 
hole we have derived the inequality (4.109). Moreover, on the horizon such a par- 


(4.143) 


ticle has been shown to have the angular velocity Qy w.r.t. a stationary observer 
at infinity, independent of the value of /. Hence with 6m = € and 6J = £ we 
obtain 
oA J 
Thus we can indeed conclude that by such adiabatic processes, adding infinitesi- 
mal amounts of energy and angular momentum at a time, the area of the black 
hole never decreases. 
Casting the relation (4.143) in the form 


Sm = = SAn + QyôJ, (4.145) 
T 


where « is given by: 


4 
k= —— Vm P, (4.146) 


mA 
it bears a remarkable resemblance to the laws of thermodynamics which states 
that changes of state described by equilibrium thermodynamics are such that for 
isolated systems variations in energy U, entropy S and angular momentum J are 
related by 


SU = TSS + QSJ, (4.147) 


where T' is the temperature and Q the total angular momentum, and moreover 
that in such processes one always has 6S > 0. The important step needed in a 
full identification of these relations is to associate the area Ay with the entropy 
and the quantity «, called the surface gravity, with the temperature: 


5S = cAn, T= (4.148) 


Sre 
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where c is an unknown constant of proportionality. Now if black holes have a 
finite temperature, they should emit thermal radiation. Using a semi-classical 
approximation, it was shown by Hawking [38], that quantum fluctuations near 
the horizon of a black hole indeed lead to the emission of particles with a thermal 
spectrum. This is true irrespective of the existence of an ergosphere, therefore it 
covers the case of vanishing angular momentum as well. The temperature of the 
radiation Hawking found was 


Ty = —. 4.149 
H 2r ( ) 


This fixes the unknown constant above to the value 


c=. (4.150) 


It seems therefore, that the thermodynamical description of black holes is more 
than a mere analogy. The issue of interpreting the entropy 

ott, fan) 
where So is an unknown additive constant, in statistical terms by counting mi- 
crostates is hotly debated and has recently received a lot of impetus from devel- 
opments in superstring theory. 

The extension of these results to the case of charged black holes is straight- 
forward. The relation between horizon area Ay and angular velocity Qy remains 
the as in (4.141), but in terms of the macroscopic observables (m, J, e) associated 
with the black hole the area is now 


Au _ 42 È g yap (4.152) 
87 2 ` ` 


The changes in these observable quantities by absorbing (possibly charged) test 
particles are then related by 


Sm = = An + QyóJ + Ode, (4.153) 
T 
where the electric potential at the horizon is 
Anep 
by = ZI (4.154) 
An 


and the constant « is now defined as 


4 
k = a Vmi — J2 — me. (4.155) 


mAH 


Again, by sending test particles into the black hole the area can never decrease: 
dAy > 0. In this way the thermodynamical treatment is generalized to include 
processes in which the electric charge changes. 
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We conclude this discussion with establishing the gravitational interpretation of 
the constant x, associated with the temperature in the thermodynamical descrip- 
tion. This quantity is known as the surface gravity, because in the spherically 
symmetric case it represents the (outward) radial acceleration necessary to keep a 
test particle hovering on the horizon without actually falling into the black hole. 

For the case of the Schwarzschild black hole this is easily established. Define 
the covariant acceleration as 


n” = Dut 
Dr 

As the four velocity of a test particle has negative unit norm: u? = —1, it follows 

that in the four-dimensional sense the acceleration is orthogonal to the velocity: 


(4.156) 


niu = 0. (4.157) 


Also, a particle at rest in the frame of an asymptotic inertial observer at infinity, 
has a velocity four-vector of the form 


1 


Combining these two equations and taking into account the spherical symmetry, 
a particle at rest must have a purely radial acceleration given by 


- 1 m 
"= —— rju = 5- 4.159 
n 2 Git =) ( ) 
Remarkably, we recover Newtons law! It follws in particular that on the horizon 
rol (4.160) 
n = p TA . 


This quantity is actually a proper scalar invariant, as can be seen from writing it 


k= In v-e | (4.161) 


rH 


as 


with ņ the invariant magnitude of the four acceleration: 


N = wn’. (4.162) 


For the case of the spinning black holes the proper generalization of this result 
gives the surface gravity as in eq.(4.155). 


Chapter 5 


Topological invariants and 
self-duality 


5.1 Topology and topological invariants 


Gravity is a gauge theory, and as such it has many similarities to Yang-Mills the- 
ories. In particular the nature of the gravitational vacuum is complicated by the 
possibilities of non-trivial topology, quantum tunneling and spontaneous symme- 
try breaking. In this section we briefly introduce some concepts and notions of 
topology. 

In the geometrical context the global topology of a manifold describes proper- 
ties characterising its connectedness, like the presence of non-contractable loops, 
spheres, etc. Clearly the way that various parts of a manifold are connected 
determines whether it can be covered by a single globally defined coordinate sys- 
tem, or requires a specific minimal number of coordinate charts glued together 
in a smooth way, i.e. such that in the overlap region of two charts the trans- 
formation from one coordinate system to the other is differentiable. Generically 
more than one coordinate chart is necessary, as in the case of a sphere in which 
conventionally the northern and southern hemisphere are mapped separately on 
the plane, with the maps overlapping in the region of the equator. Thus gen- 
eral manifolds are to be considered as glued together from a number of smoothly 
deformed regions of flat space, often in a non-trivial way. 

This has important consequences. For example, if one computes the integral of 
a certain function over the manifold, this integral must necessarily be decomposed 
into a sum of integrals over each coordinate patch. Even in the case in which 
the manifold as a whole has no boundary, the various coordinate patches do have 
boundaries. Thus, in the application of the generalized Gauss-Stokes theorem for 
the integral of a total divergence of a vector field A’,(a) over a manifold M = 
X; M;, with M; the pieces of the manifold which can be covered with a single 
coordinate chart, the integral becomes a sum of integrals over the boundaries of 
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the coordinate patches: 


V-K = f V- KÖ = f KO, 5.1 
Ja ` 2 M; ` 2 IM; tn (5-1) 


where KO (z) are the components of the vector field in the i-th coordinate 
patch, and de subscript n denotes the component normal to the boundary 0M,. 

Now if the complete manifold M has no boundary, all the boundaries 0M; are 
internal boundaries and these are integrated over twice with opposite orientation. 
Therefore, if the integrand is the same in each case, the contributions of the 
integrals over these internal boundaries cancel exactly. However, in the change 
of coordinates from one region M; to the next the value of the integrand AK 
can sometimes change and the cancellation may be spoiled. Then the non-trivial 
topology manifests itself in an apparent violation of the generalized Gauss-Stokes 
theorem: integrals over a total divergence do not have to vanish in the absence of 
an external boundary. Of course, if A, is a smooth vector field, like the physically 
observable electric field, it can not jump on the boundary 0M;: a violation of 
the Gauss-Stokes theorem requires that the vector field cannot be smooth. Only 
particular kinds of physically relevant vector fields can satisfy this requirement 
without leading to physical or mathematical inconsistencies: vector fields which 
change in a controllable way, with the change between two coordinate patches 
not leading to observable or measurable consequences. Basically such vector fields 
must be gauge fields, as these can change by unobservable gauge transformations 
on the boundary of coordinate regions. 


Example 

We illustrate these general and rather abstract arguments by a simple and well- 
known example of vector fields on the two-sphere. The standard coordinate 
system for a two-sphere of radius r is the polar coordinate system 0 < 0 < 7, 
0 < wy < 2. In these coordinates the line element on the sphere reads 


ds? = r? (do? + sin? 6 de’) . (5.2) 


From this expression it can be seen that the metric is singular —it has a zero mode 
in the y-direction— at the north and south poles, 9 = (0,7). This coordinate 
singularity is avoided by covering the upper and lower hemispheres H+, corre- 
sponding to0 < @ < 7/2 and 7/2 < 0 < 7, separately with a cartesian coordinate 
system (x+, y+), obtained for example by projection on the disc x? + y? < r°: 


z4} = rtan 5 cos P, 
for H4, (5.3) 
Y} = rtan 5 sin 9, 


and 
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z- = rcot 3 COS Ps 


for H_. (5.4) 


y- —r cot zine, 


The minus sign in the expression for y_ is necessary to have (x_, y_) form a right- 
handed coordinate system on the exterior of the southern hemisphere; it could 
be absorbed in a redefinition y —> —y. Thus the orientation of the cartesian 
coordinates is preserved in region of overlap of the two coordinate systems, which 
includes the equator 0 = 7/2: 


ty = £ =r cosg, 
(5.5) 
Y+ = —y- =rsing, 
but which actually can be extended to the entire sphere with the exclusion of the 
north and south poles by continuing (x+, y+) outside the disc x? + y? < r°. 


In the cartesian coordinates (v4, y+) the line element becomes 


4r? 
ds? = FA AF Tr E (dz? + dy?) . (5.6) 


This expression is manifestly non-singular, except in the limit (£+, y4) —> œ 
corresponding to one of the poles in each case. Observe, that the north pole is 
now a regular point of the coordinate system (x4, y+) and the south pole of the 


system (x_, y_). 
Now consider a vector field on the sphere defined in components by 


N 
Ag, = ————— (- 
t= 2a (—yt, v4), (5.7) 


where N is a normalization constant. This vector field can be represented as a 
one-form 


N 


Ay = Pope OY dvs + redyx) 
(5.8) 
N 
= 3 (+1 — cos 9) dy. 


The last expression in terms of the polar coordinates shows, that in the region of 


overlap the vector fields Az differ by an exact one-form: 


A, — A_ = Ndarctan (=) = Ndarctan (=) = Ndo. (5.9) 
x 
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It follows, that the two-form F = dA, = dA_ is unique and well-defined on the 
entire sphere: 


F = 5 Fu de” \ dx’ 


5.10) 
INr? N. ( 
= Goga py ey = 3 sinp d A dọ. 
The dual of the two-form F is defined to be the zero-form 
~ 1 N 
F = eh’ uv =, (5.11) 
2/9 2r2 


with the tensor density <“” representing the two-dimensional permutation symbol. 
As this is a constant, it is closed: dF = 0, or in components: 


OF, = 0. (5.12) 


This shows that we have constructed a non-trivial solution of the free Maxwell 
equations on the sphere, representing a constant magnetic field of strength B = 
N/2r?. This solution can be extended trivially to a static spherically symmetric 
solution of Maxwell’s equations in three dimensions, by taking the radial com- 
ponent of the vector field to vanish: A, = 0. In three dimensions this solution 
represents a magnetic monopole of strength N/2 located in the origin. 

Note that we can write the gauge transformation connecting the vector fields 
A+ in the region of overlap in terms of the elements g(y) = expi Nọ of a U(1) 
field g(x,y) in the standard form 


_ I 
A =g Ag+ 79 ldg. (5.13) 


For this U (1) field to be well-defined in the whole region of overlap of the domains 
of A+ —the whole sphere minus the poles— it must return to the same value 
after a rotation y > y+ 2r. This is the case only if N is an integer. 


The number N has a direct topological interpretation: it represents the wind- 
ing number of the map from the equator’ of the sphere to the U(1) group man- 
ifold, which is the unit circle in the complex plane. Therefore we see that the 
interpretation of the vector field Ay as a connection for a U(1) bundle (a mag- 
netic vector potential) on the sphere requires the winding number of the gauge 
transformation (the magnetic charge) to be quantized, a well-known result first 


obtained by Dirac. 


‘Note that the equator cannot be contracted to a point on the sphere whilst staying com- 
pletely inside the region of overlap where the gauge transformation is defined, as in this con- 
traction it must always pass one of the poles. 
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The example of the vector field A+ on the two-sphere illustrates how the 
naive Stokes’ theorem can be violated: even though F is a locally exact two- 


form (F = dA, in H4) well-defined and smooth over the entire sphere, and even 
though the sphere has no boundary, the integral of F over the sphere is non-zero. 
This is because F is not an exact two-form globally, and the vector field makes a 
jump on the boundary between the regions H4: 


~ 1 z 
F = J Pesa? f Pez Py 


_ (5.14) 


f F4+f F f (4; = A-) = 27N. 
Hy H C 


In the line integral the loop C = OH, = —OH_ denotes the equator, with the 
minus sign in front of OH_ deriving from the reversal of the orientation of y. 
This induces the minus sign in front of A_ in the loop integral. 

Finally we observe that a smooth deformation of the sphere will not change 
the value of the integral (5.14): if M denotes any closed surface which encloses 
the magnetic charge, and Q is the volume between this surface and the two-sphere 


S?, then 


f r- F F 
M S2 aR 


f aF =o. 
QR 


by the second Maxwell equation in the absence of magnetic charge (the Bianchi 
identity) inside Q. Geometrically we can understand this result, because the first 
line of eq.(5.14) shows that in components the integral is manifestly independent 


(5.15) 


of the metric. Hence a local change of the metric, which is equivalent to a local 
deformation, does not affect the value of the integral. Such a quantity is called a 
topological invariant. 


5.2 Topological invariants in gravity 


In gravity topology can play a role at various levels. At the macroscopic classical 
level one can consider multiply connected universes and wormholes, whilst at the 
microscopic Planck scale space-time topology may be subject to quantum fluctu- 
ations (‘space-time foam’); in analogy with other field theories like sigma models 
and Yang-Mills theories, it is expected that the quantum tunneling processes be- 
tween different topologies are dominated by finite-action solutions of Euclidean 
gravity, the gravitational instantons. 

One way to characterise topologically non-trivial solutions of the gravitational 
field equations is by the value of topologically invariant integrals over certain 
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polynomials of the curvature tensor. In four dimensions there are essentially two 
independent topological invariants of this type: the first Pontrjagin number and 
the Euler number. To define these we first introduce the dual and double dual 
curvatures with components: 


Here e“”?? is the permutation symbol, a tensor density requiring a factor ,\/g to 
guarantee the transformation character of the dual curvature as a tensor under 
local coordinate transformations. The definition (5.16) is adapted to manifolds 
with euclidean signature, but it can easily be generalized to include the lorentzian 
case. 

The first Pontrjagin number of a compact manifold M without boundary 
involves the dual curvature and is defined by the expression 


1 
pı(M) = 64r? 


whilst the Euler number is defined in a similar way using the double dual of the 
Riemann tensor: 


I, dba e RR (5.17) 


por? 


1 VR o TU 
x(M) = 555 f teva ` Sporu Rp”? Rp”. (5.18) 


For a general proof that both expressions are topological invariants it is more 


convenient to use the formulation of differential geometry in terms of the vier- 
bein and spin-connection, rather than the metric and the Riemann-Christoffel 
connection: 


R? = eeg Rivne = A, — yw, — lwp o]. (5.19) 
In this formulation, the traces in the definition of the Pontrjagin and Euler num- 


ber are to be taken over the tangent space (i.e. local SO(4)) indices (ab), and we 
can write both topological invariants in the form 


1 
327? 
where depending on the topological invariant condidered the linear operator A 
acting on anti-symmetric SO(4) tensors is either the identity or the SO(4) duality 
operator: 


cM) = I, dix 2”? Tr (Ru, ARo) (5.20) 


A RH” — R 


1 
ab, or A RH” — ze R! 


py? 


(5.21) 


respectively. From the expression (5.20) it is clear, that the square root of the 
metric has canceled inside the integral: the integral depends on the we. the 


93 


topological invariance of c(M) is now straightforward: one proves that the value 
of the integral is the same for any two spin connections which can be smoothly 
deformed into each other. Let w denote the original spin connection, and w, the 
deformed one. The difference n = wı — w is a covariant 5O(4) tensor. We define a 
path in the space of connections linking w and w, parametrized by a real number 
O<A<1: 


w(A) = w + An. (5.22) 


We show that the integral expression (5.20) is independent of A: 


de(M) 1 dR yi (À) 
= — | dev Tr | A AR (A) ) = 0. 5.23 
dÀ 167? Ja ve ( dÀ peA) (5.23) 


The last result follows, because the integrand is a total derivative which is smooth 
and well-defined everywhere if 7 is; to check this, consider the explicit form of 
the deformed curvature tensor R(X); 


Rv (A) = Re + (Dum — Dum) — y has Mul” - (5.24) 


Here the covariant derivativer of 7 is constructed in terms of the original unde- 

formed spin connection. Its derivative with respect to À is then 
d ab ab 
TEA) = Dun — Dona — 2A [nu n] 

The last term involving the commutator of two n’s is precisely the one necessary 

to complete the covariant derivatives of ņ to covariant derivatives DO) w.r.t. the 


(5.25) 


full deformed connection w(A): 


d 
L REO) = DO ny = Dy (5.26) 
Inserting this in the integral in eq.(5.23), one obtains 
de( M) — 1 4 AR ww (A) 
d 16r2 Ja orn ( dy Foo) 
1 vpo 
= -7a I, dta e” Tr [DO nA Ryo (A)| (5.27) 


1 
= = I, dix ð, (eH? Tr [ny MByo(A)]) - 


To obtain the last line we have used the Bianchi identity for the SO(4)-curvature: 


eter DR(A) = 0, (5.28) 
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which guarantees that the SO(4)-covariant derivative can be taken outside the 
trace in the integrand. As both 7 and the curvature are covariant tensors, the 
invariant 5O(4)-traces do not have a discontinuity across the boundaries between 
the different coordinate patches, and the integrand in eq.(5.27) is a globally de- 
fined total derivative on M. If the manifold has no boundary, the integral then 
vanishes, thus proving that the original quantity c(M) is independent of the 
connection indeed. 

It is of interest to note, that the Pontrjagin density can actually be written 
locally as a total covariant divergence: 


ol 


NZ du(VgJp), (5.29) 


where the current J5(x) is defined by the Chern-Simons form 


1 2 
vg Jp TEn? E”? Tr (T.a,r- — 3 PET, ) 


(5.30) 

= ao Tho opw -2 oww ) 

167? P g O 
To interpret the middle expression correctly, one should regard the Christoffel 
symbols as a set of four (4 x 4) matrices ([’,,),’, and the Riemann tensor (or its 
dual) as a set of six (4 x 4) matrices (R,,),. The trace is then taken over a 
contravariant upper and a covariant lower space-time index. In the spin connec- 
tion the trace is over the local SO(4) tangent-space indices. This result implies 
directly, that the only non-trivial contributions to the Pontrjagin number come 


from the boundaries of the coordinate patches: 


1 & , , N POO 
pı(M) = Y | dix a, ere Tr (ra,re — srororo)]. (5.31) 


16r? & Jr, BY Pe 


where the x and I denote the coordinates and connections in the i-th co- 
ordinate chart A;. That they can come from these boundaries indeed, follows 
because the current is not covariant — only the variation of its divergence is. 
Therefore we can construct another simple proof of the topological invariance of 
the Pontrjagin number: under a deformation of the metric (or the vierbein), the 
variation of the current can be written as 


S(VGIp) = pete? [Tr (RoodT,) — 8, Tr (Pp6T)]. (5.32) 
T 


or a similar form in terms of the spin connection. Because of the contraction with 
the c-tensor the last term does not contribute to the divergence of the variation 
6(,/gJp). But as the first term is covariant it does not produce a non-vanishing 
contribution to the integral by its boundary terms: its boundary terms cancel 
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pairwise between adjoining patches. Therefore the Pontrjagin number is invariant 
under smooth deformations of the manifold: 


5py(M) = 0. (5.33) 


As there is no similar Chern-Simons type current for the Euler number, this proof 
cannot be extended to that case. 


5.3 Self-dual solutions of the Einstein equations 


There is an interesting class of space-time geometries for which the two topological 
invariants — the Pontrjagin and the Euler number — are identical, or identical 
up to a sign: the (anti) self-dual geometries satisfying 


uv — 1 vo ` 4 ppv d 

R ^ = W Ros = ERM. (5.34) 
Because of the symmetry of the fully covariant Riemann tensor Rp vpo under 
interchange of the index pairs (uv) and (po), the double dual Ry,» then also 
equals the dual (up to a possible sign), and the two topological invariants are 
equivalent indeed. 

The interesting aspect of these (anti) self-dual geometries is, that they are 
automatically solutions of the Einstein equations in empty space or space-time 
—depending on the signature— with zero cosmological constant. To see this, 
construct the Ricci tensor and note that it vanishes due to the Bianchi identity: 


1 
v vÀ vpo 
R = Ra = 2g P Rurpo = 0. (5.35) 
Then automatically also the Riemann curvature scalar vanishes, and therefore 
these solutions have vanishing action, at least if the space-time manifold has 
no boundary. Euclidean solutions of this type are identified as gravitational 
instantons. 


A considerable number of such (anti) self-dual four-geometries are known, but 
only two represent regular, compact euclidean four-manifolds without boundary: 
the four-torus 74, which has a flat metric, but is non-simply connected, and the 
famous K3 manifold, which is non-flat but simply connected. Although no explicit 
expression for the K3 metric in terms of coordinates is known, its topological 
invariants can be calculated using cohomology arguments and index theory, which 
gives x( K3) = pi( K3) = 24. 

All other self-dual solutions of the Einstein equations violate at least one of the 
other assumptions: the manifold is not compact or it has a non-zero boundary. 
The earlier proof of the topological invariance of the Pontrjagin and Euler number 
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then fails: in general there can appear additional contributions from infinity or 
? In this class of solutions one finds the Eguchi-Hanson 
metric [43], and the multi-center solutions of Gibbons and Hawking [44]. 


from the boundaries. 


Taub-NUT 

A relatively simple and interesting example of these solutions is the (anti) self- 
dual Taub-NUT metric. It is the simplest example of the class of Gibbons- 
Hawking solutions. As it also has some interesting applications outside general 
relativity, I have chosen to discuss it here in some detail. 

The Taub-NUT manifold is a four-dimensional Riemannian manifold, involv- 
ing on a single scale parameter m, with Euclidean signature and (anti) self-dual 
Riemann curvature, depending on the sign of m. A convenient set of coordinates 
to describe the manifold are the four-dimensional polar coordinates (r,9,w, y) 
with ranges 0 <r<w,0<80<7,0< 0 < 7,0 <y < 27, in terms of which 
the metric is given by the line element 


Am? 


CH 
(5.36) 


The parameter m, with dimension of length, can take both positive or negative 
values. Upon factoring out m? and rescaling the radial coordinate |r/m| > r, 
the expression (5.36) becomes 


2 
ds? = (1 4 =") (ar? 4 de? + r? sin? od’) + (dw + cos Ody)’. 
r 


1 2 
— ds? = (1 + =) (dr? + 1?d6? + r° sin? Ady?) + (dib + cos Ody)” , 
r 


m2 


_4 
(7 
(5.37) 


where the sign represents the sign of m. As m? has been reduced to a multiplica- 
tive constant, one can choose units of length such that m? = 1. The right hand 
side of eq.(5.37) then directly defines the metric. In the following discussion of 
the Taub-NUT geometry this choice of metric is made, but it is straightforward 
to reinstate the general m-dependence when desired. 

In the appendix of this chapter we have collected the explicit expressions 
for the metric, the components of the Riemann-Christoffel connections and the 
covariant components of the Riemann curvature tensor. From these expressions 
one may check, that the curvature is self-dual for positive, and anti self-dual for 
negative m = —1. 

The inverse metric is encoded in the hamiltonian 


? Extensions of the construction of topological invariants including boundary contributions 
can be found in the literature, see for example the review [42]. 
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2 
Pr p 1 
2 + r(r + 2) 4 
1+- _ 
r 


which describes the geodesic flow. The canonical momenta appearing in this 
expression are derived from (5.37): 


2 2\ ; 
Pr = (1 E =) r, Pe = r? (1 E =) 0, 
r r 


4 

Py = 2 
(1 x =) 

r 


where as before an upper dot denotes a derivative w.r.t. to the affine parameter 


À. 


s+ cosb¢ = cos 0 + r°sin? 6 (144) ġ 
(ù P) Pe Py $, 


First integrals of motion can be constructed with the help of the constants of 
geodesic flow, as explained in chapter 2. Indeed, the Taub-NUT manifold admits 
a set of three independent Killing vectors which, together with the conservation 
of the hamiltonian (5.38), allow complete solution of the geodesic equations. It 
is actually convenient to express the Killing vectors as a set of four dependent 
ones, instead of three independent ones, as this leads to simplifications both in 
computations and physical interpretations. The contravariant components €©)", 
s = 1,2,3,4, of these Killing vectors can be read off from the constants of geodesic 


flow J9 = Eep, with 


J® = — sin Y pọ + DY ay — cot 8 cos Y py, J) = Py, 
sin 6 
(5.40) 


J) = cos 9 po +Z py — cot Ø sin y Py, JO = py. 
sin 


Solving these equations we get for the angular components of the four-momentum: 


po = -siny JO + cosy JO, 


pe = JO, py = J®, (541) 


whilst between the constants of geodesic flow there is the additional relation: 


J = cosð JO + sind (cos p J® + sin p J®) . (5.42) 


The radial momentum can be computed from the hamiltonian: 
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1 2 
~ 2 t [ling J® JO? 4 — (J@ — cos 9 JO ) l 
r(r +2) ((siny ane + sin? 0 ( o | (5.43) 


Before discussing some features of the solutions for the geodesics, we digress 


briefly to study the infinitesimal transformations genereated by the constants 
of motion; this gives useful general insight into the geometry of the Taub-NUT 
manifold [45]. From the last two equations (5.41) it is clear, that J@ and J@) 
generate shifts in the angles y and w, i.e. mutually commuting rotations around 
two orthogonal axes. As J?) are independent of Y, they commute with the 
rotations in the w directions generated by J“, However, they do not commute 
with J®) or with each other. A direct computation shows that their Poisson 
brackets realize the commutations relations of ordinary angular momentum: 


LJO, JO} = etk JH, i,j,k = (1,2,3). (5.44) 


As the generators J change the angle ọ in a non-trivial way, it follows that they 
do not just realize three-dimensional rotations on the submanifold spanned by the 
coordinates (0, Y). Rather, three-dimensional rotations on the spherical surface 
defined the (0, y) coordinate directions are generated by a set of related dynamical 
quantities, the restricted angular momentum three-vector I with components 


pH) = J SSP 74) 7) = JO, 
sin 6 í 
, (5.45) 
72) — JO Y ja 
sin 6 í 
The components of the restricted angular momentum three vector satisfy the 
same Poisson bracket relations (5.44) as the components of the total angular 


momentum J: 


S10, 10) = 0, (5.46) 


but with two crucial differences: (i) they induce rotations only in the (0, y)- 
directions, whilst commuting with the angle w; (ii) unlike the components of J 
they are not conserved on geodesics. This discussion of symmetries shows, that it 
is convenient to think of the Taub-NUT manifold locally as a (3+1) dimensional 
Euclidean geometry, with a three-dimensional space spanned by the spherical 
coordinates (r, 9, p) and an internal compact one-dimensional space characterized 
by the coordinate 4, but with non-trivial interplay between the flow of geodesics 
in these directions. 
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The identity (5.42) shows, that in this three-dimensional sense J“ is the 
projection of J in the radial direction, hence orthogonal to the spherical surface 
spanned by (6, y): 


(5.47) 


The contribution to the total angular momentum orthogonal to ê can be identit- 
fied as the orbital angular momentum: 


J=L 4+ JM%, (5.48) 


in the sense that 


L=(i¢ “yr x (5.49) 


Thus r- L = 0 indeed. However, in components 


LY = JO — sind cosy py = —siny po — cot 0 cos y (po — cos 8 py), 
L® = J@ — sinb siny py = cos yp» — cot 0 sin y (pe — cos 0 py), 
LE) = J®) — cosO py = py — cosh py, 
(5.50) 


and therefore the orbital angular momentum is not identical with the restricted 
angular momentum I. In fact the Poisson brackets of L do not obey the angular 
momentum algebra. Also L is not conserved, although its magnitude |L] is, 
as guaranteed by the conservation of J and J, and the orthogonality to J 
of the contributions from J“ and L in eq.(5.48). This clarifies the non-trivial 
dynamical relation between motion in the (@,y) and ¢-directions. 

Finally, as J and J are conserved separately along geodesics, it follows from 
(5.47) that the three-dimensional position vector r moves at a fixed angle w.r.t. 
the axis of total angular momentum, and the geodesic is restricted to lie on the 
cone traced out by rays from the origin around J at a fixed opening angle 


B = arccos —. (5.51) 


Obviously, this is a generalization of the motion being confined to the equatorial 
plane, corresponding to B = 7/2 (J@ = 0), in three-dimensional central force 
fields, like in the Kepler problem. Again this testifies as to the non-trivial role of 
the y-momentum py = J“) in the determination of the orbits. 
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5.4 Runge-Lenz vector and Yano tensors 


It is possible to obtain more information about the nature of geodesics by observ- 
ing that there exists a triplet of constants of motion A, i = 1,2,3, which are 
quadratic in the momenta, and which transform among themselves as a vector 
K under three-dimensional rotations [45]. This conserved vector is very similar 
to the Runge-Lenz vector in the Kepler problem. It is conserved along geodesics, 
and for J Æ 0 it satisfies the identity 


1 2 4)2 


As (K, H, J, J®) are all constant, it follows that the three-dimensional position 
vector r is constrained to lie in a plane. We have already seen, that for J 4 0 
this vector takes values on a cone with opening angle 8 (5.51) around the axis of 
angular momentum. Combining these two results it follows, that the geodesics 
are actually conic sections obtained by intersecting the cone (5.51) with the plane 
(5.52). 

It remains to construct the vector K and prove the relation (5.52). First it 
should be pointed out, that as they are quadratic functions of the momenta: 


lo. 
KO = Z KO" yp, (5.53) 


their existence is implied by a triplet of symmetric second rank Killing tensors 
KË as in eq.(2.64): 


HvV 
Du KŻ, = 0. (5.54) 


Such a triplet of Killing tensors can indeed be found. Instead of directly giving the 
explicit expressions, we construct them in terms of a set of simpler objects, known 
as Yano tensors |47]. A Yano tensor is an anti-symmetric tensor fu, = — fuu 
satisfying the Killing-like equation 


Dufur + Dy fur = 0. (5.55) 


As a result, the covariant derivative of a Yano tensor is completely anti-symmetric: 


1 
Dafa = 3 (Dafa + Dr fra + Dafo) = Dru for. (5.56) 
3 


Given one or more Yano tensors fo one can construct symmetric Killing tensors 
simply by symmetrized multiplication: 


(a l a a 
KEP = 5 (FD + FOLD). (5.57) 
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These symmetric tensors satisfy eq.(5.54) by virtue of the Yano condition (5.55). 
It turns out that the Taub-NUT manifold admits four such Yano tensors [46], 
transforming as a singlet Y and as a triplet f; under three-dimensional rotations. 
Their explicit expressions are given by the two-forms 


Y 


4(d + cos 0 dp) A dr + 2r(r +1)(r £2) sin dé A dy, 


(5.58) 


2 
fi = —6ijn 14 —) dx; A dx, + 4(dyY + cos dọ) A dz;. 
r 


A straightforward calculation shows, that the covariant exterior derivative of 
these two forms gives: 


VY = DiaYir dg” Ada” Adx* = r(r +2) sin dr A dé A dy, 


(5.59) 
Vi; = 0. 
By studying the individual terms, one establishes the stronger results 
Da Ya = DY a = DiYou, (5.60) 
and 
Da fiva = 9. (5.61) 


Thus both two-forms indeed obey the Yano condition (5.55), with the triplet of 
tensors f; actually being covariantly constant. It follows immediately that the 
triplet of symmetric tensors, defined by 


f 1 
K; w = 9 (Yò fiw + fi Yoav) ’ (5.62) 


form a vector of Killing tensors. For this triplet of Killing tensors eq.(5.52) can 
be proven. First note, that as a two-form 


r-f = —r'(r+2) sinf dð Ade + 4r(dy + cos 0 de) A dr. (5.63) 


Performing matrix multiplication and contracting with the momenta to get an 
expression of the form (5.53), it follows that 
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rK = r (fY) =r (Y f) = E (r K)” pupe 
oOo o > L 2) p? + ESIE 1 ( 0 py)? 
= Ff + -Irr pu — COS 
2(r + 2) Pr Tg Py r£2) PT ging Pe Py 
= rH 4+ (J? — J?) 
5 . 
(5.64) 
Using eq.(5.47) this gives the desired result (5.52). Note, that this equation 
remains true also for J“ = 0, in which case the geodesic flow in the three- 


dimensional (r, 0, p) space is a along a hyperbola in the plane perpendicular to 
J. If Jis along the z-axis (0 = 7/2), then the vanishing of py implies that at the 
same time wy is constant. 

The constants of geodesic flow constructed above may be interpreted as an 
analogue of the Runge-Lenz vector for the Taub-NUT geometry, because they 
can be written in the manifest three-vector form 


I 2 1 1 
K = - (1 + =) Px I Fê (a - 770?) JOJ. (5.65) 
3 UE 7 7 


Finally we comment upon the existence of the covariantly constant two-forms f;. 


They have the property that as anti-symmetric two-tensors each of them squares 
to minus the identity. In fact, they realize the quaternion algebra 


Fifi = — i; + Eijk fk (5.66) 


Thus the two-forms f; define a triplet of almost complex structures, and the four- 
dimensional real geometry of the Taub-NUT manifold can be recast in terms of 
analytic complex or quaternionic geometry. Indeed, the Taub-NUT manifold be- 
longs to the class of complex manifolds known as hyper-Kahler manifolds. An 
important consequence of being Kahler or hyper-Kahler is, that these manifolds 
have N = 2, resp. N = 4 supersymmetric extensions [48, 49, 50]. Such su- 
persymmetries give much additional insight into the properties of manifolds like 


Taub-NUT [51]. 


Appendix 


1 


03 


In this appendix I have collected the components of the connection and curvature 
tensor of the Taub-NUT manifold in the coordinate system in which the metric 
takes the form (5.37), with m? = 1. The non-zero components of the Riemann- 
Christoffel connection Pa are: 


+r] +1 
pD = , Log” — r(r ) 
rr r(r n 2) (r n 2) 
+1 2 
8 8 £ 
Peg = Te" = (r42) Top = Ty” = cot O ( 7 (r+ y) l 
r1 —2 
r.f =r %= — Tr,” S= T, 
ro er r(r +2)’ ob vo sin O(r +2)?’ 
2 
4 P cos 0 Y Y 1. cos’ 0 2 
Lro = Tor ~ (r2) Pog T e T Sny sin ð l-I? 
rT “=r p +1 r -T p 2 cot ð 
ry oyr — r(rd 2)’ ob 8 (r £2)?’ 
and 
4 1 
Loy =F (r => cos? 0 — E sin” 0, P = — sin ĝ cos 0 ( — 
> ,  __4rcos@ 9 9 2sin0 
Poy = Puy =F (r+ 2)? Poy = Puy = (r+ 2)?” 
4r 
T.” == . 
oe (r + 2)8 


ture tensor 


a 


UVR 


R 


OD, — OT sine — [T.., Ti]. . 


(r+ 


(5.68) 


From these expressions one can compute the components of the Riemann curva- 


To allow easy check of the (anti) self-duality (5.34), we give here the completely 
covariant components Rv; the non-vanishing ones are 
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Rogpoy = 


(r E 2) 


Ar? cos 0 


_2r? sin? 8 


(: 


2 cot? 0 


(r 


F 2) 


). 


in? 0 8 cot? 0 
Rrorg = = = ( 7 Soir) > 


(r2) (r +2)? 
Bos = Fy a 
Rryep = ; - =a 
Ras = — ene, 
Bovey = E 
Roveop = aiit 


(5.69) 


All other components not related by permutation symmetry of the indices are 


Zero. 
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